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We study the dynamics of quantum fluctuations which take place at the earliest stage of high- 
energy processes and the conditions under which the data from e-p deep-inelastic scattering may 
serve as an input for computing the initial data for heavy-ion collisions at high energies. Our method 
is essentially based on the space-time picture of these seemingly different phenomena. We prove that 
00 ■ the ultra-violet renormalization of the virtual loops does not bring any scale into the problem. The 

On ' scale appears only in connection with the coUinear cut-off in the evolution equations and is defined 

On . by the physical properties of the final state. In heavy-ion collisions the basic screening effect is due 

to the mass of the collective modes (plasmons) in the dense non-equilibrium quark-gluon system, 
5h I which is estimated. We avoid the standard parton phenomenology and suggest a dedicated class of 

. evolution equations which describe the dynamics of quantum fluctuations in heavy ion collisions. 
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" . I. INTRODUCTION 
> 

, In this study, we wish to approach the theory of high energy cohisions of heavy ions based on a minimal number 

g of p„.dpl.. instead of „i„g the =labo„fe .=ch..c« of facto„.a.io„-ba.ed .c.tton.g fhoo.y 0. 0„. „l...na« 

00 goal is to provide a description of the transition process which converts two initial-state composite systems (the stable 
nuclei in the "normal" nonperturbative vacumn) into the quark-gluon plasma (QGP), which is a dense system of 

O I quarks and gluons with a perturbative vacuum as its ground state. This scenario was suggested by Shuryak 20 years 

O ^' ago and it is our basic assumption that this phenomenon does indeed take place. Our main result is that this 



dense system can be formed only in a single quantum transition. This transition is similar to the process described 
by the evolution equations of deep- inelastic scattering (DIS) and differs from it only by the scale parameter inherent 
to the final state. The limit of resolution for the emission process in DIS is connected with the minimal mass of a 
jet, while in the case of the QGP final-state, the scale is much smaller and defined by the screening properties of the 
QGP itself. 

To pose this essentially quantum-mechanical problem properly, one (ideally) needs an exact definition of two main 
elements: the initial state of the system and the observables in the expected final state. 

We know very little about the initial state, and thus must rely on the following phenomcnological input only: The 
nuclei are stable bound states of QCD and therefore, their configuration is dominated by the stationary quark and 
gluon fields. The nuclei are well shaped objects; the uncertainty of their boundaries does not exceed the typical Yukawa 
interaction range. In the laboratory frame, both nuclei are Lorentz contracted to a longitudinal size i?o/7 ~ 0.1/m. 
The tail of the Yukawa potential is contracted in the same proportion. The world lines of the nuclei are the two 
opposite generatrices of the light-cone that has its vertex at the interaction point. Therefore, no interaction between 
the nuclei is possible before they overlap geometrically. 

The final state is defined more accurately. We believe that single-particle distributions of quarks and gluons at 
some early moment after the nuclei have intersected, describe it sufficiently. Thus, we may rely on a reasonably 
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well-defined quantum observable. The corresponding operator should count the number of final-state particles defined 
as the excitations above the perturbative vacuum. To develop the theory for this transition process we have to cope 
with a binding feature that the "final" state has to be defined at a finite time. This may look disturbing for readers 
well versed in scattering theory. 

Is it possible to compute the quantum-mechanical average of the operator that counts the final-state particles with- 
out precise knowledge of the wave function of the initial state? We prove that this is indeed possible, provided we have 
the data from a much simpler process, like e-p DIS. Moreover, we argue that this is possible even without constructing 
any model for the proton or colliding nuclei like, e.g., the parton model. Such a possibility is provided solely by the 
inclusive character of the measurements in both cases. Indeed, the measurement of one-particle distributions is as 
inclusive as the measurement of the distribution of the final-state electron in DIS. 

The key to our proof and to the suggested algorithm is the principle of causality in the quantum-mechanical 
measurement which has two major aspects. First, any statement concerning the time ordering must be in manifest 
agreement with the light cone boundaries. This is always guaranteed in the theory based on the relativistic wave 
equations. Second, the process of measurement physically interrupts the evolution of a quantum system, and any 
dynamical information about the quantum-mechanical evolution reveals itself only after the wave function is collapsed. 
In order to implement these principles in a practical design of a theory, one has to start with a space-time description 
of the measurement, i.e., the Heisenberg picture of quantum mechanics. 

To give a fiavor of how the method works practically, let us start with a qualitative description of the inclusive e-p 
DIS measurement (for now, at the tree level without discussion of the effects of interference). In this experiment, 
the only observable is the number of electrons with a given momentum in the final state. Something in the past 
has to create the electromagnetic field that deflects the electron. Before this field is created, the electromagnetic 
current, which is the source of this field, has to be formed. Since the momentum transfer in the process is very high, 
the current has to be sufficiently localized. This localization requires, in its turn, that the electric charges which 
carry this current must be dynamically decoupled from the bulk of the proton before the scattering field is created 
(to prevent a recoil to the other parts of the proton which could spread the emission domain). Such a dynamical 
decoupling of a quark requires a proper rearrangement of the gluonic component of the proton with the creation 
of short-wave components of a giuon field. By causality, corresponding gluonic fiuctuation must happen before the 
current has decoupled, etc. Thus we arrive at the picture of the sequential-dynamical fiuctuations which create an 
electromagnetic field probed by the electron. It is very important that the picture of the fiuctuation is dynamical. 
The short-wave Fourier components of the quark and gluon fields are more typical of free propagation with high 
momentum, than for a smooth static design of a stable proton. The lifetimes of these fluctuations are very short and 
they all coherently add up to form a stable proton, unless the interaction of measurement breaks the proper balance 
of phases. This intervention freezes some instantaneous picture of the fluctuations, but with wrong "initial velocities" 
which results in a new wave function, and collapse of the old one. An example of such a fluctuation is depicted in 
Fig. 1^, where the arrows point to the later moments of time. 

This qualitative picture has been described many times and with many variations in the literature, starting with 
the pioneering lecture by Gribov and including a recent textbook however, the sequential temporal ordering of 
the fluctuations has never been a key issue. We derive this ordering as a consequence of the Heisenberg equations of 
motion for the observables. The ordering appears to be universal only at the Born level. In general, the interference 
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effects break the sequential temporal ordering, unless the measurement prohibits the interference itself. We eventually 
arrive at a time-ordered picture which resembles the picture used by Lepage and Brodsky The difference is that 
they applied time ordering exclusively in connection with the light-cone Hamiltonian dynamics, while our picture 
holds in any dynamics as long as the expansion in terms of two-point correlators makes sense. Qualitative arguments 
supporting the strong time ordering in fast hadrons were put forward earlier by Gribov For our practical goals, 
such an ordering is significant for several reasons. 




FIG. 1. An example of fluctuation resolved by the inclusive measurement in DIS. 

First, it makes the whole picture of the evolution causal and clearly indicates that the inclusively measured fluctu- 
ations are independent of the nature of the inclusive probe. 

Second, the causality allows one to prove important dispersion relations. These relations provide a tool for the 
renormalization which is also subjected to causality, viz., there cannot be radiative corrections to a stable proton (or 
nucleus) configuration. This condition is physical and fully equivalent to the on-mass-shell renormalization of the 
asymptotic states in a standard scattering theory. Thus, we do not treat the intermediate quark and gluon fields as 
fundamental fields that have certain states. Only hadrons (in DIS) or collective modes (in dense matter) can be the 
states of emission. 

Third, the causality in the inclusive measurement requires that the quark and gluon fields, which might exist after 
the measurement, are fully developed as the virtual fluctuations before the measurement. This is true both for freely- 
propagating hard quarks and gluons and for their final-state hadronic equivalents (jets). Therefore, the fiuctuations 
have to satisfy the condition of emission. This condition requires that (i) the field of emission exists as a detectable 
state, and (ii) that two different states can be resolved. Unlike QED, QCD has an intrinsic criterion of resolution 
which requires that the states of emission in the physical vacuum must be hadrons. Thus, e.g., there cannot be a 
quark or gluon emission which carries away an invariant mass tuq less than several hadronic masses. In other words, 
the modes with wave-length A > jtt-q ^ ^ lGev~^ cannot propagate. The process of formation of the final-state jets 
has to be fully accomplished in a few fm. This general rule, however, must be changed if a quark or gluon becomes 
a part of a dense system, like the QGP, where the boundary of the soft (and hopefully still perturbative) dynamics 
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is determined by a parameter similar to the plasmon mass (or the Debye screening length) and not a hadronic scale 
. This condition selects fluctuations that are materialized as the QGP after the coherence of the nuclei is broken. 
Fourth, our analysis, based on the causal picture of the evolution, indicates that the only scale that may appear 
in the theory is connected to the physical properties of the final states. This scale penetrates into the theory only 
in connection with the infra-red (coUinear) problem but not in connection with ultra-violet renormalization. This 
observation allows one to outline the scenario of formation of the dense matter at the first Ifm of a nuclear collision 
as follows: The nuclei probe each other locally at all scales which are consistent with the spectrum of the allowed 
final states. The most localized interactions create the final-state fields with the highest pt which are the least 
sensitive to the collective properties of the final state. (A similar process takes place in DIS at sufficiently high 
3> "^"hadr ^%CD-) Thcsc statcs are formed as (quasi-) particles by the time r ^ l/pt- In this domain, the 
structure functions and their rates of evolution must be very similar for both nuclear collisions and DIS. Therefore, 
the DIS data obtained with precise control of the momentum transfer may serve as an input for the calculation of the 
nuclei collision. The states with lower pt are formed at later times. This docs not much affect the QCD evolution in 
DIS, where the phase-space of the final states is almost unpopulated and only the shape of the jets may lead to the 
power corrections. In the case of nuclei collisions this is not true, as the low-pj states are formed in a space which is 
already populated with the earlier formed (virtually decoupled) fast particles . These fast particles produce screening 
effects. The process of the formation of the new states saturates at pt ^ mjj, where mo is an effective mass of partons 
in the (highly non-equilibrium) collective system formed by the time tjj ^ 1/mjj. To quantify this scenario one has to 
derive evolution equations that respect specific properties of the final states in the QGP. A draft of such a derivation 
is given in the last two chapters of the paper. 

To conclude, the standard inclusive e-p DIS delivers information about quantum fluctuations which may dynamically 
develop in the proton. We do not need to know the "proton's wave function" (since the proton is a true fundamental 
mode of QCD, this notion is not even well defined). The theory can be built on two premises: causality, and the 
condition of emission. The latter is also known as the principle of cluster decomposition which must hold in any 
reasonable field theory. What the "resolved clusters" are is a very delicate question. These states should be defined 
with an explicit reference as to how they are detected. Conventional detectors deal with hadrons and allow one to 
hypothesize about jets. QGP turns out to be a kind of collective detector for quarks and giuons. 

A vast search for a scale in ultra-relativistic nuclear collisions has been undertaken by McLerran et al. This study 
started with the formulation of the model for a large nucleus as a system of valence quarks moving at the speed of 
light (McLerran- Venugopalan model Q). The original idea of this approach was that the scale is associated with the 
density of valence quarks, that is, with the parameters of the initial state of the nuclear collision. The model was 
gradually improved (see ||^, and references therein) by the accounting for the radiative corrections with emphasis on 
the domain of small xp and a derivation of the BFKL equation, which is known to have no scale parameter. 

The idea that screening effects should be taken into account, even at the earliest moments of a collision of two 
nuclei has been articulated earlier and with different motivations by Shuryak and Xiong Q and by Eskola, MuUer 
and X.-N. Wang . In Ref. 1^, the cut-off mass (expressed via the temperature) was introduced to screen the 
singularities of the cross-section of the purely gluonic process, gg ^ ng, in the chemically non-equilibrated plasma, 
g'^rigiuon/T . In Ref. ||^, the non-equilibrium Debye screening was introduced in order to prevent the singular 
infrared behavior of the hard part of the factorized cross-section of the mini-jet production. The expression for the 
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screening mass in B is very similar to our equation (10.16). However, in all of these papers, the origin of the structure 



functions of the colliding nuclei is not discussed, and there still remains an uncertainty connected with the choice 
of the factorization scale. We go one step further and prove that the screening parameters of the collective system 
determine the coUinear cut-off self-consistently; these screening parameters are already in the evolution equations for 
the structure functions of individual nucleons (and therefore, in the nuclear structure functions). 

II. OUTLINE OF CALCULATIONS, MAIN RESULTS, AND CONCLUSION 

The QCD evolution equations are the main focus of this study. Historically, their development started along two 
qualitatively different lines which eventually merged into what is currently known as the DGLAP equations. Despite 
their common name, there are still certain concerns whether the marriage is happy. The one-to-one correspondence 
between the DGL jl^ and the AP |ll[] equations is evident only at the Born level (without the virtual radiative 
corrections). A discussion of what happens in the higher orders and particularly, what the argument of the running 
coupling is, still continues (see, e.g. Refs. |p^ , |l3| ). The variations appear due to ab initio different approaches. While 
the AP equations re-interpret the equations of the operator product expansion (OPE) based on the factorization and 
the renormalization group (RG) for e-p DIS, the DGL equations are derived as a probabilistic picture of the multiple 
process. The AP equations just accept the running coupling of the RG for the exclusive Feynman amplitudes, while in 
the DGL scheme one has to calculate the radiative corrections as constituents of a real process explicitly and assemble 
them into an effective quantity (the running coupling). 

The first claim about a closer correspondence between the DGL and the AP equations, at least in one loop beyond 
the Born approximation, was made in Ref. |l^ (widely known as DDT) almost 20 years ago. The ultimate goal was to 
recover the running coupling in the natural environment of the QCD evolution ladder by collecting various radiative 
corrections into an effective form-factor instead of relying (as does the OPE method) on the RG for the Feynman 
amplitudes of the S-matrix scattering theory. This program was a success only for one subgroup of ladder diagrams, 
vis., when the gluon is an offspring of the fermion line and the dynamics is driven solely by the splitting kernel Pqq. 
This was a smart choice, and success was almost predetermined by the well known fact that in physical gauges the 
renormalization of the coupling is fully defined by the gluon propagator. Nevertheless, the DDT study did mighty 
battle with the double-logarithmic terms which naturally appear in radiative corrections. It was a formidable task to 
show that the running coupling depends on the highest transverse momentum in the splitting point. Neither DDT, 
nor anybody else since that time, attempted to do the same for other types of splittings, though it is not obvious if 
the same result can be obtained for the coupling constant which accompanies the kernel Pqg (when a single gluon 
line carries the lowest of three transverse momenta). We intend to bridge this gap and assemble the running coupling 
constant of all other kernels. The kernel Pgg (when all three lines are gluons) is of extreme interest since one may 
expect competition between three gluon self-energies with all different momenta. The case of pure glue-dynamics 
provides an ultimate test for the understanding of the dynamics of QCD evolution. We address it first, and use it 
later as a touchstone for all the tools we design. 



For the sake of completeness, we begin in Sec. [H with the basics of the so-called Keldysh-Schwinger technique [Il5| 



in the form adjusted for the calculations of the inclusive observables [16|. The net yield of this section is summarized 



in Eq. (3.36) for the fixed-order two-point field correlators. This equation shows that the correlator T>{x,y) has two 
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contributions: from the initial data, D{^, ry), and the source, n(^, ry), acting at an earlier time. In Sec. [II C we discuss 
the connection between Eq. ( |3.36 ) and the spectral representation of field correlators and explain why the initial data 
can be dropped. In Sec. |^ these equations are used to derive the equations that lead to the picture in Fig. |^ and 
prove that the photon of e-p DIS must be treated as part of the proton. 



We also examine (in Sec. IV) what types of fluctuation are active in the e-p DIS which is inclusive with respect to a 
quark jet (without observation of the final-state electron). In this problem, the photon "belongs" to the initial-state 
electron and the measurement is affected both by QED and QCD fluctuations which may develop in the initial-state 
electron. 

In Sec. ^ we derive the evolution equations of glue-dynamics yet without the radiative corrections. Even at this 
level, we discover the second independent element of the evolution, vis., the longitudinal fields, which are absent both 
in the DGL and the AP equations. [| The reason is that, in our derivation, the evolution equation naturally appear 
in a tensor form (^.3|). Projecting this equation onto the normal modes we obtain two interconnected evolution 



equations, (5^) and (p.q), for the transverse and longitudinal (static, non-propagating) gluon modes. This should 
not be considered a surprise. Indeed, it is well known that the photon on the top of the evolution ladder (just proved 
to be a part of the proton) has transverse and longitudinal polarizations which are probed by the structure functions 
F2 and -Fl, respectively. It is not clear, a priori, why the gluon field should be only transverse. If the longitudinal 



field is eliminated from Eqs. (5.4) and (5.5) by a fiat, the remaining single equation is exactly the DGLAP evolution 
equation of glue-dynamics. Its kernel is singular, and it is common to regulate it with the aid of the (+)-prescription. 
In Ref. p8[ this prescription which incorporated self-energy corrections (in the fixed coupling approximation) was 



derived. For a similar singularity in Eq. (5.5) this method does not work. Moreover, it is not clear what kind 
of coupling accompanies the transition between the longitudinal and transverse fields in the evolution equations. 
Potential consequence of this new-type of infrared divergence might be a disaster since the standard remedies of 
scattering theory are not designed to handle any problems related to the longitudinal fields. In scattering theory the 
latter are never treated as observables. The specifics of heavy-ion collisions are diff'erent because the observables must 
be defined at a finite time of evolution. Unlike scattering theory, the longitudinal fields cannot be absorbed into the 
definition of the final states. Before looking for a cure, we must make sure that this is not a false alarm, e.g. perhaps 
the infrared divergences do not cancel between the different rungs of the emission ladder. Thus, we need to study the 
full set of radiative corrections in the order g'^. 

A full set of radiative corrections includes different types of diagrams. Two diagrams in Fig. ^ are virtual vertex 
corrections. Their vertex loops are formed by three T-ordered propagators, -DqOj in ^00 (diagram (a)) and three re- 
ordered propagators, Dn, in Vm (diagram (b)). Even though all propagators which are connected to these vertices 
are retarded, the vertex loops have a priori arbitrary ordering of their space-time arguments. This may hurt the nice 
causal picture of the lowest order approximation. Our main observation is that the net yield of the sum of these 
two diagrams does not change if the vertex Vboo is replaced by the retarded vertex Vret (diagram (c)) and the vertex 
Viii is replaced by the advanced vertex Vadv (diagram (d)). In this formulation, the picture of the evolution remains 



^It was Jianwei Qiu who first paid attention to the special role of the static part of the quark and gluon propagators 
(he called them special or contact propagators) in computing the twist-4 corrections to the DIS process. However, the main 
emphasis of this study was on reproducing the high-twist effects in the perturbative part of the OPE-expansion. 
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totally causal even in the order g^. 

The retarded properties of the new vertex functions in coordinate space can be translated into the analytic properties 
in momentum space. To comply with the light-cone dynamics and the light-cone gauge = 0, wc must treat the 
coordinate a;"*" as the temporal variable. The conjugated variable is the light-cone energy p~ . We prove that the vertex 
function Vret{p^ k — p,—k) is analytic in the upper half-plane of complex p~ . The light-cone energy p~ corresponds 
to the latest momentum resolved in the course of the evolution (which has the lowest component p"*" of the light- 



cone momentum and is, in fact, probed in the inclusive measurement). These results are obtained in Sec. VI B and 
Appendix 1. Knowing the analytical properties of the vertex function we are able to prove the dispersion relations 
which allow one to find the real part of the vertex via its imaginary part. It occurs that all singularities which may 
affect the final answer are collinear and already appear in the calculation of the imaginary part which is connected with 
real processes. Thus, these singularities are physical and do not change in the course of a calculation of the dispersion 
integrals. Besides, they are always connected with emission into the final states and never with the absorption from 
the initial states. 

The dispersive method provides a tool for the UV-renormalization of the vertex function, viz., the subtractions in 
dispersion relations. These subtractions must be made at a value p^ = —Q, oi the light-cone energy in the domain 
were the imaginary part is zero. This domain occurs to be limited to a finite segment ( |6.23| ) of the real axis in the 
complex plane of p^ . The position = — of the subtraction point can be transformed into the value of p^ , vis., 
p^ — — /i^ = —p^n — pi- Under the assumption of strong ordering of the transverse momenta, pt » kt, the allowed 
values of are < fi^ < p^. Even though this limitation is a direct consequence of causality, it may seem abnormal 
from the RG point of view. Indeed, we cannot take an asymptotically high renormalization point in p^. However, in 
a causal picture of gradually increasing resolution along the line of evolution, this limitation is absolutely natural. It 
is illegal to renormalize what is not yet resolved! 

If we proceed with renormalization at an arbitrary scale /i^, then the calculation of the vertex function leads to 
the terms log^{p^/fi^) along with log{pl/iJ^) . The double- logarithms find no counterparts among either virtual 
loops or real emission processes. Therefore, we immediately undermine the leading logarithmic approximation (LLA) 



in which the DGLAP equations are derived (see Eq. (7^)). Besides, in order to have log{p^/iJ?) >> 1 one should 
choose a very low renormalization scale /i^ which would contradict all beliefs about perturbative QCD. These problems 
compel us to look for a physically-motivated condition for the renormalization. Surprisingly, this condition is readily 
found at the opposite end of the allowed interval. If we take /i^ = p^ which corresponds to the subtraction point 
at p^ — 0, then all of the above problems miraculously disappear. The vertex function, as well as all retarded 
self-energies (also calculated via dispersion relations), become free of logarithms, log{pl/ ji^), which originate from 
ultraviolet renormalization (where they are most expected to come from). Technically, the toll is that we have to 
look for the large logarithms in real processes. However, from the physical point of view, this is a gain rather than a 
loss, since we get rid of the logarithms which are just an artifact of the renormalization prescription. In fact, we can 
promote the status of the running coupling; instead of being a formal attribute of the renormalization group it may 
become a dynamical form-factor due to the real processes. Before we search for these large logarithms, let us try to 
understand why the subtraction point = is a natural physical point for the renormalization. 

According to the principle of causality, the full picture of the QCD evolution develops before the measurement. 
Thus, the system of fields must be coherent. It represents, at most, virtual decomposition of a hadron in terms of an 
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unnatural set of modes which emerge as propagating fields only after the interaction happens. Therefore, the radiative 
corrections to the propagation of fields should not cause any phase shifts along the line of hadron propagation and 
hence, the real part of virtual loops should vanish for the state of free propagation. This is exactly what is achieved 
by the on-mass-shell renormalization of the asymptotic states in scattering theory. If the hadron is assumed to move 
along the light cone, the real parts of the vertices and self-energies should vanish at — 0. 



Contributions ( [7.9D and (7.10) of the retarded gluon self-energy and vertex (both rcnormalized at the point — 0) 
to the evolution equation are double-logarithmic with respect to the coUinear cut-off e. This cut-off already appears 
in the imaginary parts of the self-energy and vertex. Therefore, it is connected with real processes. Even if the 
divergent terms do cancel, one still has to understand why. The physical mechanism of this cancelation can only be 
due to interference, and it has to work even at the classical level. The collinear divergence is both a classical effect and 
an artifact of theoretical model with massless charged particles. Therefore, it should not be treated as an artificial 
phenomenon caused by spurious poles of a special gauge. Indeed, the physical mechanism of the collinear divergence 
can be easily translated into the language of space-time. 

Consider a charged particle that instantaneously changes its velocity from vi to V2 at a time t — 0. At some later 
time t, the picture of the field is as follows. Outside the sphere of radius R ^ t, there is an unchanged static proper 
field of the initial-state particle. Inside this sphere, there is the newly created proper field of the final-state particle. 
The field of radiation (which has the old and the new static fields as the boundary conditions) is located on the sphere 
itself. If U2 < 1 then the charge is inside the sphere and decoupled from the radiation field. At V2 = 1, the charge never 
decouples from the radiation field and this leads to the divergence in the integration over the time of the interaction. 
Thus, assigning a finite mass m to the field is a test for true collinear divergence. This is also a cure, since then V2 < I 
and the problem disappears. Alternatively, one may provide at least a single interaction with a third body which 
interrupts the emission process after some time At. In both cases, we encounter a dimensional physical parameter 
which regulates the time-length of interaction between the charge and the radiation field. If m 7^ 0, then the smallest 
of two quantities. At and m^^, is the physical cut-off for the emission process. Let us assume that some mass mo 
can be introduced as a collinear cut-off. If this cut-off cancels in the full assembly of diagrams in any given order, 
then the theory has no single dimensional parameter. This is exactly the case of pQCD where the parameter Aqcd 
establishes, at most, an absolute boundary of its validity. All practical calculations are possible only for momenta 
higher than some Qq ^ Aqcd- All the divergences, including soft and collinear, are usually removed by means of the 
dimensional regularization with reference to the KLN theorem; however, this is not the case we address. 



In Sec. VIII, we calculate contributions of real emission processes to the causal QCD evolution. The diagrams in 
Figs. ||a and^b, with the unitary cut through the vertex, and the interference diagram in Fig. do not lead to terms 
with log{pf) . They, however, include many terms with the collinear cut-off which eliminates a part of phase-space 



which does not satisfy the emission criterion. In Sec. IX we study the radiative corrections in the LT-transition mode. 
They confirm the general trend. 

Thus, we have proved that in the time-ordered picture of evolution the UV-renormalization should be performed 
in such a way that radiative corrections do not affect the initial state and bring no scale into the problem. The scale 
is brought in by the emission process and is defined by parameters of the final state that emerge after the collision. 
The most important conclusion that follows from this observation is that the evolution equations must be modified, 
viz., the universal running coupling that accompanies splitting kernels in the DGLAP equations must be exchanged 
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for the form-factors which depend on the type and scale of the emission process. This is required by causahty of the 
inclusive measurement. 

It is high time to understand why we need a theory with a physical cut-off of the collinear emission and why we 
cannot rely on the standard approach of pQCD. The latter was constructed as a theory of the S-matrix at asymptotic 
energies. The concept of asymptotic freedom is formulated precisely in this context. The running coupling is always 
introduced as a certain assembly of radiative corrections regardless of whether the perturbation series is divergent. In 
pQCD, the way to assemble vertices and self-energies is defined by the structure of Feynman amplitudes of exclusive 
processes, and it is very important that the collinear (mass) singularities are, at any price, removed from this assembly. 
The most popular technical tool to do this is based on dimensional regularization of the Feynman amplitudes. This 
method does not distinguish between UV-, IR-, and mass-singularities, despite the different physical nature of these 
divergences. Sometimes, one may trace that cancelations of certain divergences, are due to interference of various 
partial amplitudes. In other cases, one may observe that the soft amplitudes exponentiate and die out in the quantities 
like total cross-sections. The underlying calculations always implicitly rely on the structure of the space of states, 
which includes arbitrarily soft states of emission (otherwise, the exponentiation makes no sense). Such a space of 
states is completely physical in QED but not in realistic QCD. Usually, the collinear singularities are forced to cancel 
with the aid of the KLN theorem which redefines the states of scattering by admixing a probabilistic distribution of 
soft (collinear) states. The proof of the KLN theorem is intimately connected with the definition of the states, vis., 
it is sensitive to the parameter of resolution, which is hidden in the definition of the degenerate states. Once again, 
the theorem was derived under the influence of QED where even very soft photons can be resolved in the asymptotic 
states, provided a sufficiently long time is allowed for the measurement. To approach this idealized regime in QCD, 
one has to limit the theoretical analysis to exclusive processes with very energetic final-state jets (insensitive to several 
additional soft hadrons which might simulate the emission of extra soft gluons). However, even in QED there are cases 
when the recipe of KLN is not supposed to work. For instance, in the QED plasma, the space of states is populated, 
and we are not free to assign arbitrary statistical weights to the collinear states. 

In nuclear collisions, the space of states has a quite different structure. Unlike a scattering problem, the states must 
be defined at a finite time of evolution. Furthermore, the longitudinal fields cannot be absorbed into the definition 
of stable charged particles (by means of on-mass-shell renormalization) and arc not geometrically separated from the 
fields of emission. Nevertheless, the technical problems are not frightening, since the normal modes in dense and exited 
(hot) nuclear matter (to the first approximation) can be computed perturbatively. In Sec. we compute the mass of 
a transverse plasmon (a collective mode), which is now the final state for the emission in the QCD evolution process. 
The distribution of gluons that induce a plasmon mass is defined by the process of multiple-gluon production, and the 
plasmon mass is the function of the time of the distribution measurement. The result of this calculation clarifies the 



correspondence between the two approaches. It is summarized by Eq. (10.14): The larger the transverse momentum 
Pt of the emission is, the earlier it is formed as the radiation field, the less it is affected by the interaction with other 
partons, and the smaller its effective mass is. This is exactly the limit when the OPE-based calculations of DIS and 
the here advocated self-screened evolution should merge. The waves with a very large pt just do not "see" the medium 
and have to hadronize into jets as if there was no medium. In order to quantify this observation, one has to estimate 
the corrections to the LLA on both sides. These corrections lead to the broadening of jets in DIS and radiation losses 
of fast partons in AA-coUisions. 
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Formula (10.14) was derived in three steps. We started (in Sec. XA) with an attempt to compute the mass of 
the transverse plasmon in the nuU-planc dynamics, which we essentially relied on in the derivation of the evolution 
equations. It turns out that the dispersive part of the retarded self-energy is always proportional to and may 
contribute only to the renormalization of the propagator. At first glance, the non-dispersive part seems to be capable 
of generating the effective mass of the plasmon. However, even in the case of a finite density of emitted gluons, the 
non-dispersive part remains quasi-local because of the singular nature of longitudinal fields in the null-plane dynamics, 
and it vanishes in the course of renormalization. We conclude that it is impossible to generate an effective mass of a 
quasi-particle in the null-plane dynamics. 

In order to smooth out the geometry of the longitudinal fields, we had to choose another Hamiltonian dynamics 



which is described in Sec. XB . This new dynamics uses the proper time r, =t^ — z^, as the time variable. For each 
slice in rapidity, r is just the local time of the co-moving reference frame. An advantage of this dynamics is that it 
naturally incorporates the longitudinal velocity of a particle as one of its quantum numbers and visualizes the process 
of the particle formation by the time ^ 1/pt. In its global formulation, the proper-time dynamics uses the gauge 
= 0. The physical mechanisms of the screening effects (like the generation of the plasmon mass) are known to be 
localized in a finite space-time domain. This understanding allowed us to use the gauge = in a local fashion. In 
order to estimate the plasmon mass we approximated this gauge by the temporal-axial gauge (TAG) = of the 
local (in rapidity) reference frame. 



An estimate of the plasmon mass is accomplished in Sec. XC. Once again, in local TAG only the non-dispersive 
part of the retarded self-energy (in which the propagation is mediated by the longitudinal field) and the tadpole term 



(with the inertia- less contact interaction) lead to the plasmon mass, which is given by Eq. (10.14). The low-pt mode 
of the radiation field acquires a finite effective mass as a result of its forward scattering on the strongly localized (and 
formed earlier) particles with qt ^ pt- The smaller pt is, the larger the plasmon mass is. Thus, the whole process, 
that converts the nuclei into the dense system of (yet non-equilibrium) quarks and gluons, must saturate. In this way, 
the scenario approaches the end of the "earliest stage", i.e., the dense system is already created and the collisions 
between the partons-plasmons start to take over the dynamics of the quark-gluon system. 

To derive the evolution equations with screening, we had to find a space-time domain where the proper-time 
dynamics is close to the null-plane dynamics. The limit is found at large negative rapidity, viz., in the photon 
fragmentation region. In the last section of the paper, we show that when the effective mass of the plasmon is used as 
the pole mass of the radiation field in the evolution equations, the evolution equations become regular. The coUinear 
singularity of the splitting kernel Pgg as well as a similar singularity in the LT-transition mode are screened by the 
plasmon mass. Thus, we obtain a simple version of the evolution equations with screening, which can be applied to 
the fiuctuations in the photon fragmentation region. They become the DGLAP equations when s oo. In order to 
obtain the evolution equations valid in the central-rapidity region, we have to consider the whole problem in the scope 
of the "wedge dynamics" |2^j2^. The last thing we estimate within our "local approximation" is that the plasmon 
mass leads to strong suppression of all radiative corrections in the order aj.. 

The net yield of this paper is as follows. The interaction between the two ultrarelativistic nuclei switches on 
almost instantaneously. This interaction explores all possible quantum fluctuations which could have developed by 
the moment of the collision and freezes (as the final states) only the fiuctuations compatible with the measured 
observable. These snapshots cannot have an arbitrary structure, since the emerging configurations must be consistent 
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with all the interactions which are effective on the time-scale of the emission process. In other words, the modes of 
the radiation field which are excited in the course of the nuclear collision should be the collective excitations of the 
dense quark-gluon system. This conclusion is the result of an intensive search of the scale inherent in the process 
of a heavy-ion collision. We proved that the scale is determined only by the physical properties of the final state. 
Eventually, the scenario for the ultra-relativistic nuclear collision promises to be more perturbative than the standard 
pQCD. It will be free from the ambiguities of the standard factorization scheme inherent in the cascade models. 

This approach is novel for the field of ultrarelativistic nuclear collisions, but the underlying physics is not new. 
Quantum transitions in solids always excite collective modes, like photons in the medium, phonons, or electrons with 
effective masses and charges (polarons). Similar phenomena are known in particle and atomic physics. For example: 

1. Let an electron-positron pair be created by two photons. If the energy of the collision is large, then the electron 
and positron are created in the states of freely propagating particles and the cross-section of this process accurately 
agrees with the tree-level perturbative calculation. However, if the energy of the collision is near the threshold of the 
process, then the relative velocity of the electron and positron is small, and they are likely to form positronimn. It 
would be incredibly difficult to compute this case using scattering theory. Indeed, one has to account for the multiple 
emission of soft photons which gradually builds up the Coulomb field between the electron and positron and binds 
them into the positronium. However, the problem is easily solved if we realize that the bound state is the final state 
for the process. We can still use low-order perturbation theory to study the transition between the two photons and 
the bound state of a pair ||T^ . 

2. Let an excited atom be in a cavity with ideally conducting walls. The system is characterized by three parameters: 
the size L of the cavity, the wave-length A <C i of the emission, and the life-time At = l/F of the excited state. The 
questions are, in what case will the emitted photon bounce between the cavity walls, and when will the emission field 
be one of the normal cavity modes. The answer is very simple. If cAt <^ L, the photon will behave like a bouncing 
ball. When the line of emission is very narrow, cAt 3> L, the cavity mode will be excited. It is perfectly clear that in 
the first case, the transition current that emits the photon is localized in the atom. In the second case it is not. By the 
time of emission, the currents in the conducting walls have to rearrange charges in such a way that the emission field 
immediately satisfies the proper boundary conditions. We thus have a collective transition in an extended system. 

From a practical point of view, these two different problems are united by the method of obtaining their solutions. 
A part of the interaction (Coulomb interaction in the first case, and the interaction of radiation with the cavity walls 
in the second case) is attributed to the new "bare" Hamiltonian which is diagonalized by the wave functions of the 
final state modes. This allows one to reach the solution in a most economical way. The modes of the proper-time 
dynamics are introduced with the same goal of optimizing the solution of the nuclear collision problem by an explicit 
account of the collision geometry (adequate choice of the quantum numbers) already at the level of normal modes. 

A list of known examples can continue. The example of nuclear collisions is new in one respect only. Traditionally, 
the collective modes are constructed against the existing "heavy matter background". In nuclear collisions, the 
collective modes are created simultaneously with the matter which supports them. 
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III. THE EQUATIONS OF RELATIVISTIC QUANTUM FIELD KINETICS 



The equations of quantum field kinetics were derived in Ref. ||16[ with the goal of calculating the observables which 
cannot be reduced to the form of the matrix elements of the S- matrix, viz., the composite operators that cannot be 
written as the T-ordered products of the fundamental fields. The inclusive rates (or inclusive cross-sections) which 
are really the number-of-particles operators, are quantities of this sort. In we concentrated on the phenomena 
associated with the dynamics of fermions. In this study, the main emphasis is on the detail calculations in the gluon 
sector, and thus, we use equations of QFK for vector gauge field. 



A. Basic definitions 



The calculation of observables such as inclusive cross sections is based on work by Keldysh ||T5[ . It incorporates a 
specific set of exact (dressed) field correlators. These correlators are products of Heisenberg operators, averaged with 
the density matrix of the initial state. For the gluon field they read 

T^T^i^, y) = ~i{A--'^{x)A^-'P{y)). y) = -*(A^^^(j;)A'^^"(a;)), 

Bfo''^{x,y) ^ -^{T{A{xr■'-A'>■'0{y))), D?J^"^(:r, y)= -*(rt(A'^^"(x)A''^'3(y))), (3.1) 

where T and are the symbols of the time and anti-time ordering. They may be rewritten in a unified form, 

^AB^ix^y) = -i{T,{A----{xA)A''--^{ys))), (3.2) 

in terms of a special ordering Tc along a contour C = Cq + Ci (the doubled time axis) with T-ordering on Co and 
T^-ordering on Ci. The operators labeled by '1' are T^-ordered, and stand before the T-ordered operators labeled by 
'0'. Recalling that 

A(x) = S''T{A{x)S) = T''{A{x)S^)S, (3.3) 



we may introduce the formal operator Sc — 5*^5, and rewrite (3.3) using the operators of the in-interaction picture, 

Bfj,"^{x,y) = -t{TM^--^{xA)A'-^''{yB)S,)) • (3.4) 

Except for the matrix form, the Schwinger-Dyson equations for the Heisenberg correlators remain the same as in any 
other technique. An elegant and universal way to derive them (that does not rely on the initial diagram expansion) 
can be found in Ref. . For the gluon field these equations are of the form 

Dab = Dab + J2 ° ° + J2 ° ° ^sb, (3.5) 

R RS 

where the dot stands for both convolution in coordinate space, and the usual product in momentum space (providing 
the system can be treated as homogeneous in space and time). Indeed, the only tool used to derive these equations 
was the Wick theorem for the ordered products of the operators. (The type of ordering is not essential for the proof 
of the Wick theorem [^). The second term on the right corresponds to the interaction with the "external" field A. 



The formal solution to this matrix equation (3.5) can be cast in the following symbolic form. 
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[D-%b = [D-%b-Uab, (3.6) 



Explicit expressions for the self-energies Hab emerge automatically in the course of the derivation of the Schwinger- 
Dyson equations; we find, 

1 



R,S=0 

-^(-1)^+^ E (-1)^+" /rf€rf^V^f4,,,,(^,x,y,^)D^-^ . (3.7) 

i?,s=o 

where the first and the second term correspond to the loop and tadpole diagrams, respectively. The 3-gluon vertex 
(in coordinate representation) is defined as 

v,t/%p(-,y,-) - (-i)^+^+^ ^P"(;;;^)]^^^' . (3.8) 

g 6Aiy{zp) 

The coordinate expression for the bare 3-gluon vertex is 

(^) ^ (^1, ^2, X,) = -gfabcb'^^id'l - d]) + gP-'id^ - d^) + 5^"(af - dl)] , (3.9) 

where one has to put xi — X2 — = x aX the end of the calculation. In the momentum representation, it may be 
written as 

ylBC;abc{Pl^P2,P^) = ~W 5 ab5 AC r'^^^Hpi - P2V + g^^ {P2 - P^T + g'^ (p^ ' Plf] ■ (3.10) 



The next (third) order correction to the three-giuon vertex is obtained by means of the Eq. (3.7), after differentiating 
the second-order self-energy with respect to the field A, 

^ bcf.RSB\^^y^^) — \ ^) <:Aff \ ~ 

g SAHzb) 

= ~igV{xR)DRs{x, y)Viys)DsB{y, z)VizB)DBRiz, x) , (3.11) 

where all Lorentz and color indices in the second equation are dropped for brevity. 

In this paper, the tadpole term in the gluon self-energy will not be considered beyond the one- loop approximation. 
Hence, we need only the bare four-gluon vertex, 

ylsRS^abrs = - 9^ SabSarS AS [fnar fnbsig"^ g"" " S""/") + 
+ fnasfnbr{g"''g"' " ff"""/") + /nafc/„.. (ff""/" " 9"" 9^')] ■ (3.12) 



The four types of operator ordering which enter Eqs. (3.1) are not linearly independent, i.e., there exists a set of 
relations between the field correlators and the self-energies: 

Doo + Dn=Dio + Doi, Hoo + Hn = -Hio - Hoi . (3.13) 

These indicate that only three elements of the 2x2 matrices D, 11, etc. are independent. To remove the over- 
determination let us introduce new functions 

Dret = Dqo — Dqi, Dadv — Dqq — -DiQ, Di = I?oo + Du; 

n,.et = noo + noi, n„d„ = Hoo + nio, ni = noo + nii, (3.14) 
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One possible way to exclude the extraneous quantities is to use the following unitary transformation [IE 

1 



D = R-^DR, n = R-^UR, R 



V2 



1 1 

-1 1 



In this new representation, the matrices of the field correlators and self-energies have a triangle form, 



D = 






^adv 


n = 


Hi 




Dret 




^adv 






(3.15) 



(3.16) 



Applying the transformation ( 3.16| ) to the matrix Schwinger-Dyson equations (3.5) we may rewrite them in the 
following form: 



"^adv — Do^dy ~\~ Dadv ^ n-adi; ^ '^adv i 
Dl ~ Di + Dret O Ilret O Di + Di O Ilarf,, O Dadv + Dret O Hi O Dadv ■ 

Using the definitions (3/7), ( |3.14 ) and the identities, 



(3.17) 
(3.18) 
(3.19) 



Don{xi,X2)DoQ{yi,y2) - £'oi(a;i, 2^2) £'10(2/1, ^2) = -Dio(a;i, a:2)£'oi(yi, 2/2) - £'ii(a;i, a:2)£'ii(2/i, 2/2) = 

= ^[-Dret(a;i,a;2)£'i(2/i,2/2) + £^i(a;i, a;2)£'Qrf„(2/i, 2/2)] 



(3.20) 



we may write the explicit expressions for the retarded and advanced gluon self-energies in the approximation of the 
bare vertex, 



n 



af3,ab / 



i^^y) = -iTC^x)D;';jfax,2/)Kr'(2/)D( ''^'{y,x) 



4^"/- ^-^^^Cf,) 



c'bf 



+Vp{x)Df'^-'{x,y)Vf,;,{y)D'^J^\y,x)] + 6{x - y)Ii'i!^^\x,x) 



(3.21) 



It is easy to see that one of the propagators in the loop keeps track on the temporal ordering of the arguments of 
these correlators. The second propagator, is nothing but the density of physical states with the on-mass-shell 
momentum which mediates the propagation of the gluon. 




FIG. 2. The retarded gluon self-energy. The arrows point to the latest temporal argument. The crossed lines correspond to 
correlator Di which represents the density of states of free propagation. 
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In order to find the tadpole part of the retarded gluon self-energy, Iltdpi{x, x), we have to go back to the Schwinger- 
Dyson equation which contains the following fragment, 



The sum over the contour index R involves two Green functions, Doo{z, z) and _Dii(z, z), with coinciding arguments 
which are ill-defined since the arguments are assumed to be ordered. Only the sum Do(){z, z) + Dii{z, z) — Di(z, z) 
is unambiguous. Using the following chain of transformations (which is very similar to ( 3.20| )), we obtain 



(3.22) 



Since the correlator Di in the last equation is on mass shell, and D^et on the left is not, the second term on the right 
drops out and the tadpole term of the retarded self-energy becomes 



n 



ap.ab 



(3.23) 



tdpl V" 1 ~ / 2 

This term is entirely due to the instantaneous contact interaction of the propagating field with on-mass-shell modes. 



B. The formal solution of the integral equations 

The solution we shall look for now is equivalent to the rearrangement of the perturbation series for observables. 
This rearrangement is useful when specific features of the non-equilibrium system must be taken into account. First, 
let us introduce two new correlators: 

Do = Dret - Dadv = Dw - Dqi , (3.24) 
which coincides with a commutator of the gluon fields and thus disappears outside the light cone, and 

Ho = Ilret — ^adv = ~ Hqi , (3.25) 

which is the commutator of two charged currents and has the same causal properties as ( |3.24 ). The integral equation 



for Do may be derived by subtracting Eq. ( p.lTl) from Eq. (|3.18| ), 

Do = i^O + Dret O Ilret O Dq + i^O ° ^adv O ^adv + Dret O XIq O Dadv 



(3.26) 



The sum and the difference of Eqs. (3.19) and ( 3.26| ) give corresponding equations for the off-diagonal correlators Diq 
and I?oij 



Doi = Dai + Dret ° ^ret ° Doi.lO + Dm O Uadv ° J^adv - Dret ° Hoi O Dq, 



10 10 



dv • 



Since Eq. (3.17) for the retarded propagator may be identically rewritten in the same form, 

Dret = Dret + Dret ° ^adv ° 'Dadv + Dret ° ^ret ° D,.et — Dret ° ^adv ° T^adv, 



(3.27) 



(3.28) 



we may use Eq. (3.27) and derive the corresponding equations for the T- and T^-ordered propagators. 
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Do 



Doa + Dret ° ^ret O Doo + Doo O Uadv O T^adv + Dret O IIii O Da, 



dv ■ 



(3.29) 



This chain of routine transformations reduces the equations that make up all elements of the matrix correlator Dab 
to a unified form. On the one hand, this representation shows that linear relations between the correlators (or, 
explicitly, different types of orderings) hold even for the equations that the correlators obey. On the other hand, this 
representation of the equations singles out the role of retarded and advanced propagators over all other correlators. 



In order to understand why their role is special, let us transform them further, and begin by rewriting of Eq. (3. IS) 
for the density of states Di identically as 



(1 - Dret O Hret) O Di = Di O (1 + Uadv ° ^adv) + Dret O Hi O T>adv ■ 

Since T>ret ° ^ret ° Dret — ^ret ~ Dret, it is casy to show that 

(1 + Dret O nret)(l - -Dret O Ilret) = 1 • 

Furthermore, we have the following two relations: 

(1 + ^adv ° J^adv) = D^Q^ O Dadv , (1 + D^et O Uret) = D^et O U^^J , 



(3.30) 



(3.31) 



(3.32) 



where D^^J^^x) and I?^p^(x) are the left and the right differential operators of the wave equation, respectively. They 
explicitly depend on the type of the Hamiltonian dynamics, including the gauge condition imposed on the field A{x). 



Multiplying Eq. ( 3.30 ) by (1 + D^et ° ^ret) from the left, we find the final form of the equation. 



Di = Dret O £'(0) ° ^1 ° ^(0) ° ^<^dv + D^et O Hi O Dadv 



(3.33) 



Repeating these transformations for the other equations, we obtain the corresponding forms that are most convenient 
for the subsequent analysis: 



Dio = 'Dret O D, . O Dio O O Dadv - 'Dret ° IIio O Dadv, 

01 01 01 



(3.34) 



Do 



Dret O Dff^] O Doa O D,^^ O Dadv + Dret ° Hii O Dadi;- 



'(0) " ^11 " "^(0) " ■^"'f ^ '^ret ^ iiii ^ ^adv (3.35) 

While equations (3.17) and ( 3.18| ) are the standard integral equations which have an unknown function on both sides; 



however, after transformation the unknown function appears only on the l.h.s., and we may consider Eqs. (3.33)^(3.3?;) 
to be the formal representation of the required solution. 

At this point, the first and the most naive idea is to ignore the arrows indicating the direction the differential 
operators act in, and to rewrite Eqs. (3.34)-( 3.35| ) in the momentum representation. Then the first term in each of 



Eqs. ( 3.34 ) will contain the expression like p'^S{p^), which equals zero. This feature reflects the simple fact that the 
off-diagonal correlators Diq and Dqi are solutions of the linearized homogeneous Yang- Mills equations. However, this 
approach does not appear to be sufficiently consistent: we lose the identity between Eqs. ( 3.34 ) and ( |3.26 ), and make 



it impossible to generate the standard perturbation expansion in powers of the coupling constant. Eqs. ( 3.35 ) will be 
corrupted as well. 

A more careful examination of Eqs. ( |3.34 )-(3.35) shows that all four dressed correlators Dab can be found as the 
formal solution of the retarded Cauchy problem, with bare field correlators being the initial data and the self-energies 
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being the sources. Indeed, integrating the first term of each these equations twice by parts, we find for all four 
elements of Dab to be 

01 J 01 J 01 

Doo(x,y) =JdU[pcin';;>^I)ret{x,Od''DooitvW^adviv,y)~ y"dV^/Dret(x,C)[±i^o"' +nii(e?7)]D,rf„(77,y) . (3.37) 



The equations (3.36) and (3.37) are the basic equations of rclativistic quantum field kinetics. They are identical 
to the initial set of Schwinger-Dyson equations, but have the advantage that the time direction is stated explicitly. 
There are two terms of different origin that contribute to any correlator (and, consequently, to any observable). The 
first term retains some memory of the initial data. The length of time for which this memory is kept depends on the 
retarded and advanced propagators. The second term describes the current dynamics of the system. A comparison 
of these two contributions allows one to judge if the system has various time scales. 

C. Initial data and spectral densities. 

We perform all calculations in the gauge — which is complementary to the null-plane dynamics with the the 
time variable . In this gauge, the on-shcU correlators and the propagators of the perturbation theory are as follows, 

d^^ip)^-,^''^?^^^^. (3.38) 

ad. P+{P ±10)-Pt P+ 

The polarization tensor 11^'' appears only between the retarded and advanced propagators, i.e., in the combination 
[Dret{p)^{p)Dadv{p)]'^^ ■ Both propagators contain the same projector, d^^''{p) which (by the gauge condition) is 
orthogonal to the 4- vector . So, of the general tensor form, only two terms survive, 

lV^''{p)^g^"' P^ wi{p)+P^p''w2{p) + ... . (3.39) 

The others, like p^n'^ + n^p'^ or n^'n"^ , cancel out. The invariants wi and W2 can be found from the two contractions, 

- d,,.(p)n^"(p) = 2/ wi(j}), and n^n,n''''(p) = {p+f W2{p) , (3.40) 

independently of the other invariants accompanying the missing tensor structures. (The projector df^,j{p) is defined 



below in Eq. (3.42).) The solutions of the Schwinger-Dyson equations (3.17) and (3.18) for the retarded and advanced 



gluon propagators. 



can be cast in the form. 



>->(r.t AP) = 77^^ + —+72 7fl^ • (3-41) 

p-^+p^wY'{p) > l-w^^^'ip) 

The first term is the propagator of the transverse field and has a normal pole corresponding to the causal propagation. 
It contains the polarization sum for the two transverse modes, 
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cT (p) = ~d^^P{p)d;{p) - -g^^'' + ^ -^-^ - P'j^, (3.42) 



which is orthogonal to both vectors n'^ and p^. The second term in Eq.(3.41) is the response function of the longitudinal 
field which does not lead to causal propagation. When radiative corrections wi and W2 are dropped, the bare 
propagators are recovered with a clear separation of the transverse and longitudinal parts. This separation tells us 
that the longitudinal field in the null-plane dynamics has only the component A" which is generated by the charge 
density of the source current which enters Gauss' law but not the equations of motion. If all other components of the 
current vanish, then (by the current conservation) we have d+j^ = and, consequently, j+ — j^{x^) = j^{x'^ — x^). 
Therefore, the purely static field of the null-plane dynamics is produced by the charge traveling at the speed of light 
in the positive x^-direction. The radiation field is everything else minus this static field. Physically, the field is called 
the radiation field if it geometrically decouples from the source. Therefore, we can formulate a criterion for radiation 
in the null-plane dynamics. The radiation field has to be slowed down with respect to the static source. There is 
no problem with real massive sources, they are never static in the afore-mentioned sense. However, if the source is 
put into the infinite momentum frame, a paradox arises and the problem re-appears through the poles in the 

propagators. Fortunately, in the case of QCD the radiation is always massive, which resolves the paradox and screens 
the light-cone singularity in the emission process. However, since the screening is a non-local effect associated with 
the longitudinal fields, it cannot be derived in the singular geometry of the null-plane dynamics (see Sec.^). 

/R\ ^ /R\ fR\ /R\ 

With the shorthand notation, W^"^' (p) — p^ + p^ wl^' {p), and (p) = 1 ^ ^2 {p)i we readily obtain. 



which also exhibits a clear separation between the transverse and static fields. The first term on the right side 
is extracted by contraction with the metric tensor g'^'^; the second one is extracted with the aid of the projector 
p^p'^ / [p'^Y ■ The combination ( |3.43| ) enters the right side of Eq. (3.34) and represents the density of states into which 
the giuon field may decay. 



The first term in Eq. (3.34) corresponds to the initial data for the field correlator. The role of this term is two-fold. 
According to the picture qualitatively described in the Introduction, the final states in the inclusive measurement are 
fully developed being the virtual fiuctuations before the measurement. Therefore, we encounter two types of initial 
fields. The fields of the first type (with the correlators labeled as D*) are real quark and gluon fields of the proton or 
nuclei before the collision. They always have sources and no initial data which can be expressed in terms of the free 



fields. Hence, the first term in Eq. (|3.34 ) for these fields has to be dropped. To understand further consequences of 



this step, let us assume that the full retarded and advanced propagators in Eq. (3.34) are renormalized according to 



some condition (to be specified later). At the renormalization point, the retarded and advanced self-energies vanish 



and, therefore, HretDfJ. = D.^/Dadv — 1 ■ This leads to 



'(0) - ^(0) 



Dto = ZDl^- D;,,H^oD:,„, (3.44) 



where Z is the residue of the pole of the renormalized propagator. Thus, the absence of the initial data for the 
free gluon field is equivalent to the requirement that its renormalization factor Z equals zero. Hence, this is not a 
fundamental field and its spectral density is formed solely by multiparticlc states into which it decays. (This option 
to make a distinction between the elementary and composite fields was intensively discussed by Weinberg [p6[.) 
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The fields of the second type (with the correlators labeled as D*) are used to decompose the colliding composite 
objects in terms of the true normal modes of the final state which are excited only after the collision. For example, in the 
case of the two nuclei colliding, the normal modes are plasmons which have effective masses m. To some approximation 
(connected to the possible damping of excitations in the collective systems) we have, Z?]^ = D"^ ^ S{p'^ — to^)6'(p°). 



and we may transform Eq. (3.34) into the known expression for the spectral density, 

/>oo />oo 

T>*'"'{p^)= dm^p^"'{m^)S{p^ ~^m^) = Z d'"'{p)S{p^ -m^)e{p°)+ dm^cr^"' {m^)S{p^ ^ m^)9{p°) . (3.45) 

^0 "'0 

Masses that enter this spectral representation shield the abundant coUincar singularities that may appear in the 



evolution equations. 



IV. TEMPORAL ORDER IN INCLUSIVE PROCESSES 



The program for computing the quark and gluon distributions in heavy-ion collisions relies on the data obtained in 
seemingly more simple processes, e.g., ep-DIS. It turns out that differently triggered sets of data may carry significantly 
different information. In this section, we discuss two examples and show that even minor change in the way the data 
are taken may strongly affect what is actually observed. 

Let us consider a collision process where the parameters of only one final state particle explicitly measured. Let 
this particle be an electron with momentum k' and spin a' . Deep inelastic electron- proton scattering is an example of 
such an experiment. All vectors of final states which are accepted into the data ensemble are of the form a^, (k')|X) 
where the vectors \X) form a complete set. The initial state consists of the electron with momentum k and spin a 
and the proton carrying quantum numbers P. Thus, the initial state vector is aJ.(k)|P). The inclusive transition 
amplitude reads as (X|a(j'(k') S aJ.(k)|P) and the inclusive momentum distribution of the final-state electron is the 
sum of the squared moduli of these amplitudes over the full set of the non-controlled states \X). This yields the 
following formula, 

^ = (P\aAk)S^alik')a^,ik')Saiik)\P) , (4.1) 

which is just an average of the Heisenberg operator of the number of the final-state electrons over the initial state. 
Since the state \P) contains no electrons, one may commute electron creation and annihilation operators with the 
^-matrix and its conjugate S't pulling the Fock operators a and to the right and to the left, respectively. Let 
ip^\x) be the one-particle wave function of the electron. This procedure results in 

^^ly f dxdx'dydy'^l'JixWitlix'm^^ ( 1^ ) (y)^!^, (y') . (4.2) 

dk' 2f^J y y v^k^ V ;vk<T V n ^s-i'{x)d^{y)\6^{y')S-^ix')J^ /^ka \y)^k a J v ; 

Since the electron couples only to the electromagnetic field, then, to the lowest order, 

^ = ^ E/ dxdyi^itUy)-^i'-Jix){P\ ^^-<\x) ^^^\y)\P)^'C\y)^\^l{^) > (4-3) 

where A^''^ [x) is the Heisenberg operator of the electromagnetic field. Already at this very early stage of calculations, 
the answer has a very clear physical interpretation. Since only the final-state electron is measured, the probability 
of the electron scattering is entirely defined by the electromagnetic field produced by the rest of the system evolved 
from its initial state until the moment of interaction with the electron. 
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The correlator of the two electromagnetic fields in Eq. (|4.3|) is the function iD^Q''(x, j/). Therefore, we may use 



Eq. dpI) to express it via the electromagnetic polarization tensor of the proton. Summation over the spins of the 



electrons brings in the leptonic tensor, L^y(fc,A;'). If q = k — k' is the space-like momentum transfer, then the DIS 
cross-section is given by 

where ttJ'q (q) is the notation for the correlator of the two electromagnetic currents, 

<;(x,2/) = (P|r(x)r(y)|P) . (4.5) 
(Alternatively, we may obtain the answer by means of the Yang-Feldman equation |p3| , 

A(x) = j AD(l|(a;,z/)j(2/) , 

where j(?;) is the Heisenberg operator of the electromagnetic current and D^g)(x,?/) is the retarded propagator of 
the photon.) The correlator of the currents is the source of the field which has scattered the electron. Here, both 
photon propagators are retarded and indicate the causal order of the process. The measurement analyzes all those 
fiuctuation in the hadron, which have developed before the scattering of the electron and created the virtual photon 
probed by the electron. 

The tensor tt^q (x, y) is naturally expressed via the two quark correlators, 

<(2^,2/) = Tr7^Gio(x,2;)7''Goi(y,x) . (4.6) 

The correlators Gio and Gqi are evolved from the earlier times via equations (3.26) and (3.27) of Ref The field 
correlator generated via the source of the previous step of evolution is 

Gqi — » Gqj^ = — GretSgj^Gadu ■ (4-7) 

The newly created correlation Gio in the state of free propagation is connected with the density of the final states, 

Gio ^ Gf^ = Gff^ - GfgjEfgGf^^ . (4.8) 

Similarly, the self-energies S(5i are defined via the correlators Goi and Dio which can also be causally evolved from 
previous stages of the evolution. In this way, the causal ladder of the fluctuations is built even without the concept 
of partons. 

Early discussions of the role of the light cone distances in high energy collisions resulted in Gribov's idea of 
the two-step treatment of inelastic processes p5[ ; the gamma-quantum first decays into virtual hadrons and later 
these hadrons interact with the nuclear target. This idea looks very attractive since it explains the origin of the two 
leading jets in electro-production events, corresponding to the target and the projectile (photon) fragmentation. It 
also provides a reasonable explanation of the plateau in the rapidity distribution of the hadrons. However, this elegant 
qualitative picture contains a disturbing element, i.e., the way the words "first" and "later" are used. 

In order to understand how Gribov's process may be observed practically, let us change the observable in the same 
deep inelastic process initiated by the ep-interaction. Let us trigger events on the high-pt quark or gluon jet in the 
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final state regardless of the momentum of the electron. This means that all states that are included into the data set 
are of the form, a\{p)\X), where a\{p) is the Fock operator for the final state quark. The corresponding observable 
is the number of the final-state quarks. To the lowest order of perturbation theory, this is a process of the type 2 — > 1. 
Its inclusive probability is 

^ = (PKik) S^a{ip)a,{p) S at (k)|P) . (4.9) 
After commutation of the Fock operators with S and 5^ we arrive at the expression, 



dp 



(4M) 



With reference to the Bogolyubov's form of the micro-causality principle p^ , which (in its simplest form) reads as 

JIr-^ f^^I^) = 0' -y?>0 and x"^ > y\ (4.11) 

we may argue that the space-time points x' and y' (where the quark jet is created) are inside the forward light cone 
of the possible points x and y of the electron scattering. To the lowest order, the graph for this process is depicted at 
Fig. |. 




FIG. 3. Fluctuation which is "active" in the inclusive measurement of jet distribution in DIS. 

The photon propagator is retarded. To produce the quark, the hadronic target has to be hit either by the photon 
coming from the first step of the virtual fragmentation of the electron or by one of the partons formed by the further 
virtual fragmentation. By causality, both the electron and the hadronic target must develop appropriate fluctuations 
before the moment of the creation of the final-state quark. Thus, the change of the trigger drastically affects the 
information read out of the data. Triggering on the high-p^ jet in the same process allows one to filter out fluctuations 
corresponding to Gribov's picture. 



V. THE EVOLUTION LADDER TO THE LOWEST ORDER 



The way to derive evolution equations in the causal picture of QFK has been shown in Refs. |1^,|27|. In the case 
of pure glue-dynamics, we start with the Eq. (3.7) for the tensor Il^on, and iterate this equation with the aid of the 
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Eqs. ( ^3.34 ). To the lowest order in g^, the vertices are considered to be bare (without the loop corrections). The 
iteration process explicitly accounts for the two types of fields corresponding either to the core of the proton (these 
are not yet touched by the evolution) or to the fields of anomalous fluctuations (which eventually become the field 
of radiation.) Therefore, every off-diagonal correlator, Dqi or Dio, corresponds either to the initial field of colliding 



composite particles, or to the final state of the gluon field. In the first case, we use equations (3.34) and express the 
field correlators via their sources formed during the preceding evolution. 



a 



dv • 



(5.1) 



The first term in Eq. (3.34) is dropped, which corresponds to the main physical assumption that the proton or nucleus 
is a fundamental mode of QCD and that there is no free gluon fields defined by the initial data. In the second case, 
we replace Dio by the densities of the final states, 



Dio 



D* = D* 



D* n*D* 



(5.2) 



Here, the first term corresponds to the immediate creation of the final state on-mass shell gluons (therefore, their 
propagators have no radiative corrections). By keeping this term intact, we implicitly recognize that the states of 
freely propagating color fields do exist. This is both the main assumption and the corner-stone of perturbative QCD. 
Though this step allows one to approach the exclusive high-energy processes from the "safe" side of the QED-like 
perturbation theory, it completely disregards the existence of the hadronic scale and creates artificial problems for 
collinear processes. Further analysis indicates that in the picture with realistic final states, this term has either to 
be modified by the redefinition of the normal modes {e.g., in dense matter), or be abandoned (if the final states are 
hadronic jets). The second term corresponds to multiparticle emission (via the intermediate off-mass-shell gluon field) 
and is more physical. 

For now, we limit ourselves by the condition that only one two-point correlator of this expansion can originate from 
the QCD evolution of the initial state and that the emitted gluon is on-mass-shell, 

n^i» = -^ / -0^V^T/{P,k-P,-k)[Iiretimo^{k)Ba,.{k)^^^^^^ . (5.3) 

This is the full tensor form of the evolution equations of pure glue-dynamics to order g'^. They have to be projected 
onto the normal modes of the gluon field and be rewritten in the form of scalar equations for the invariant functions. 
Our next goal is to find the evolution equations for the invariants, wi and W2 of the gluon field. Using the Eqs. ( 3.40| ) 



and (3.43) we extract evolution equations for the various invariants from the tensor evolution equations (5.S). We 
obtain two equations for the invariants of the gluon source: 



(P) 



NcUs 



d'^kS+[{k - pY' 



W?{p)Wt{p) "'^ 'w^{k)Wt{k) 
>Vf(p)Wj4(p) W^{k)W^{k) 



(5.4) 



,,(01) 



ip) 



2Ncas 



W.^{p)Wiip) TT^ 

where Pgg{z) is a well-known splitting kernel, 



d'^kS+[{k-pf 



(1/z- 1/2)2 Pw(°^^(fc) 
W^{p)W^{p) Wf (fc)Wj4(/c) 



(5.5) 
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+ z{l-z) 



(5.6) 



z 1 — z 

Within the accuracy of this approximation, viz, when '-^■'lloi ^ is used as the basis for the iteration procedure, the 
propagators yV^'{p) and 'W^{p) on the right-hand side do not include loop corrections. 



The first term on the right side of Eq.(5.4) has a transverse source in it, which describes a step of evolution when 
the initially transverse field pattern remains transverse after the emission. In what follows, we shall refer to this type 
of transition as the TT-mode. The second term in this equation contains a longitudinal (static) source and after one 
act of real emission the new field pattern becomes transverse. This field can be, e.g., the proper field of the static 
source, which became a propagating transverse field after the static source has been accelerated (TL-transition mode) . 
The second equation describes the process of creation of the new static field of the source after the acceleration has 
been terminated (LT-transition mode) . There is no LL-transition mode since any rearrangement of a charged system 
between two different static configurations requires at least two emissions. The first emission extinguishes the old 
static field, the second emission creates a new one. 

Connections between the invariants w^^{p) and the structure functions follow from the description of measurement. 



dp 



_ ip'^p wi{p) dG{x,pj) 



dp 



ip+ W2{p) 

W^{p)W^{p) 



(5.7) 



Wf[p)W^{p) dp} ' , 'y2KFI'^2 

where c is some common normalization constant. In terms of these functions the evolution equations acquire a habitual 
form, 



dG{p+,pl) 



dp} 



dk^dkf 



_S+\{k-pfl dG{k 

dp — Pgg(z) 



p. 



dk} 



dp-'-^^^{^-lroik-,k^) 



(5.8) 



GiP^,p't) 



dk^dkf 



dp-5+[{k^pf 



1-2 dG{k+,kf) 



(5.9) 



We must account for the static components in course of the QCD evolution as long as we wish to describe this process 
as developing in time, and as long as the fields are considered not only in the far zone, but also in the near zone, 
where they continuously interact with the sources. 



None of the integrations dp in Eq. (5.S) is singular at the point A;+ = p^ . Indeed, assuming the strong ordering 
of the emission, pt^kt, we obtain, 



dG{p+,pi) _ N,as 



dp} 



dk+ 
2k+ 



1 



Pgg{z) / dkt 



dG{k+,k 



dkf 



Pt ^ 



(5.10) 



It is common to regulate the singularity of the splitting kernel Pgg[z) at z = 1 with the aid of (+)-prescription. 
In Ref. 1 18 1, this prescription to order g^ was obtained using two diagrams depicted in FigjH The singularity in the 
one-loop correction to the free propagation compensates for the divergence due to the coUinear emission. 
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FIG. 4. Two diagrams used to form the (+)-prescription. 

In the case of emission (that converts transverse field into the longitudinal one), we encounter the same type of 
end-point singularity. Integration dp^ in Eq. ( [5.9| ) is singular and leads to the pole (1 — 2)^^ at fc+ =p^, i.e., 

gi:,,p^)^-^N^J^^ —--———G{k+), (5.11) 

where 

is the x-fraction of the glue converted into radiation, integrated over all transverse momenta starting from some low 
(confinement ?) scale Qq. 



Though the pole in Eq. (5.11) is similar to the one in Pgg{z) no self-energy correction to this process exists, since 
the two-point function cannot convert a T-mode into an L-mode. The only hope for the possible removal of this 
divergence is connected to the radiative corrections of the emission process. Besides, unlike (transverse) partons the 
longitudinal fields are beyond the factorization scheme and the RG for the S-matrix. It is not at all obvious what 
kind of running coupling accompanies the splitting kernel. Therefore, we have to consider the evolution with the 
elementary cell at least to the order g"^. 



VI. THE EVOLUTION EQUATIONS TO THE ORDER 



s 



There are two kinds of problems which cannot be resolved within the approach of a simple ladder. First, we 
encounter coUinear singularities which may be (at least partially) compensated by loop corrections like W^'^{p). 
Second, the loop corrections and the additional real emission processes are required to assemble the running coupling. 



Let us start with the full expression for ^'*)noi(p), which emerges after the vertex of the g'^-order (3.11) is substituted 



into the Eq.(3.7) 



R,S=0 

xV,%'^y)I)%^^'^^''\y,z2)VZ::^{z,)Bl^^^^^^^ . (6.1) 



In this equation, we still keep all the gluon correlators dressed. We decompose it using the following guidelines. First, 
we divide terms into two major groups corresponding to R — S and R ^ S. The former gives rise to the virtual 
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vertex corrections, while the latter describes all possible real emission processes. The loops of virtual corrections are 
formed by the diagonal propagators. They are ultraviolet-divergent and require renormalization. Inside each group, 
every off-diagonal correlator, Dqi or Dio corresponds either to the initial field of the colliding composite particles, or 



to the final state of the gluon field. Thus, we again use Eqs. (5.1) and (5.2) along with the previous agreement that 
only one two-point correlator of this expansion can originate from the QCD evolution of the initial state. 

In this approximation, two types of radiative corrections appear in the evolution equations. First, we have to 
calculate the basic polarization tensor Hoi to order by accounting for the vertex corrections. Second, we should 
account for the loop corrections in all propagators. The formal counting of the order comes from the following 
decomposition, 

Doi = — DreilloiDad^ = —Dret [1 + ^retDret] Hqi [1 + Dadv^adv] Dadv 

« -Dretl^^'^Uoi + ^^^UretDret^^'^Uoi + ^^^UoiDadv^^^Iladv]Dadv , (6.2) 

where ^"^11 ~ g". After separation of the transverse and longitudinal modes, we obtain, 

p2 ^01 (p) p2 (4)^01 (p) ^ (2)^|(p)(2)^01(p) 

^ 2p , (6.3) 



' + 2"'"S(p)'^'-nP)- (6.4) 



A. The gluon self-energy corrections 

Various gluon self-energies '■^^nyis(p) of the order g"^ are needed to compute the elementary cell of the QCD ladder 
to order g^. Though the results are well-known, we consider it instructive to go through their derivation in some 
details. 

The simplest of ^^''Hab, the self-energy (^)nf^(p— fc), is the correction to the emission process due to the "mass" of 
the decoupled gluon which splits further into two massless gluons. These two gluons belong to the space of the final 
freely-propagating states. Therefore, 

n?o^(p) = J d*qe{q+)0{p+ q+)S{q+q- q^,)5[{p+ q+){p- q- ) [pt - gl)^] 

xV^7(l~P^P^-<l)dpa{q)V''''''{q,-p,p-q)d,,,{q~p) . (6.5) 

After integration over with the aid of the first 5-function we obtain the following expressions for the transverse 
and longitudinal invariants wf^^ and wf^'^: 

J d\-'tS[{p+ - q+){p~ -^)- {pt - qlf] . (6.6) 

The divergence of this integral can be viewed in different ways. First, let us notice that the integration d'^qt involves 
only a (5-function and the result of this integration is 7rj3+/g+. Thus, the divergence becomes connected with the dq^ 
integration, and requires both lower, and upper cut-offs in the remaining integral for w^^^{p), 



#10 
,,#10 



ip) 
ip) 



ig^Nc 



dq^ 



-p'P,9iiT) 



2U-^ 
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<'(p) - --^e{p')9{p+mp\^) , (6.7) 

where we have introduced the functions 

/3i(p+,e) = -27V,^' '^P^^(^) ^ Ar^(^^_41n^), when e^O, (6.8) 

P2{p'-,e)=N,J' "^(^l-2^y ^ ^, when e^O. (6.9) 

Being introduced in this way, the cut-ofF may be thought of as a regulator for the spurious poles in the gluon 
correlators Dio{q) and I?oi(9 ^p) which form the loop. However, such an interpretation is excluded by the fact that 
these correlators are on-mass-shell and include only transverse fields. Therefore, the poles of the polarization sums 
d{q) and d{q — p) are due to the transverse momenta of the on-mass-shell loop gluons and the singularity of the 
remaining integration is a physical collincar (mass) singularity. To remove this singularity, one has to eliminate the 
collinear region by a fiat, viz., by introducing a physical parameter of resolution for the emission process. The value 
of this parameter has to be extracted from the data, e.g., from the shape of jets in DIS. In the case of heavy-ion 
collisions, the cut off will naturally come from the finite density of the final-state particles. Indeed, in the coordinate 
space, a collinear singularity is caused by the infinite time-of-interaction between the radiation and its massless source. 
A single interaction with the "third body" suffices to interrupt the emission process and to allow the field of radiation 
to decouple from its source. 

If the cut-off parameter is sufficiently small, than the integration results in the well known numbers associated with 
the Gell-Mann-Low beta-function. However, one should keep in mind that if the Ine-tcrm is not miraculously canceled 
due to some interference process and has to be kept finite, then the magic "(lliVc/3)" becomes an approximate number. 
In fact, the cut-ofF e can be re-expressed via the lower limit of the invariant mass of the two final-state gluons and it 
has to be even kept in the finite quantity f32- 

To the lowest order, the retarded self-energy of a gluon can be computed via its imaginary part. This statement is 
not trivial. Indeed, the retarded self-energy is as follows, 

+yfa7il - V.V. -g)Dr' ''^(g)<;,(g, - 9)Dg„y'(g - v)\ ■ (6.10) 

Though in both terms we encounter the retarded propagators, only their transverse parts are truly causal. Thus, if 
we wish to use the dispersion relations to calculate the real part of H^et (which is responsible for the phase shifts in 
the propagation of the gluon field) via the imaginary part (which is directly connected with the real processes), we 
have to separate causal and non-causal parts in Eq. ( B.lOj ), 



D,,M)^DVJM + D^'^\^). D,M^D^aliq) + D^^\q) . (6.11) 
Thus, we obtain the dispersive part of the retarded (advanced) self-energy, 

+V{q~p,p,-q)D^{q)V{q,-p,p-q)D^^l^{q-p)] , (6.12) 
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and the non-dispersive part, 



n^Z'^hp) = I ^[V{q-p,p,-q)D^'^H<l)V{q,-p,p-q)D^{q-p) + 



i:t)^' 4 7 (27r)4' 

+V{q-p,p,-q)Di{q)V{q,-p,p-q)D^'^\q-p)] . (6.13) 

The difference, 

Ho = Ilret - liadv = H^S^ - H^^] = Hqi - Hio , (6.14) 

is the imaginary part of the retarded self-energy, which has been aheady computed. Projecting onto the two normal 
modes, we obtain the real part of the invariants as the dispersion integral of the imaginary part with one subtraction 
at some value p~ = —fl , 



Re<"*(p) = w'Jp) = - I dw- ( 

■' TV J \(jO~—p~ UJ~ - 



Imu;;:^*(a;) , (6.15) 



where 



Imw^/^Lo) = ^e{p+uj- - p2)sign(p+)/?,(p+, e) . (6.16) 

An explicit calculation of the dispersion integral results in the renormalized retarded self-energy, 

Re<*(p) = -l^Pjip+,e) In , (6.17) 

where the subtraction point p~ = — O is translated into p"^ = —^j? = — p+fl — 

The non-dispersive part identically vanishes for the real part of w"* (because of the polarization properties of the 
three-gluon vertex). For the non-renormalized non-dispersive function we obtain a divergent integral, 

which is independent of p~ and is, therefore, a quasi-local function proportional to (5(a;+). It identically vanishes after 
one subtraction at any value of p~ . Thus, after the renormalization, the retarded self-energy is totally causal both for 
the transverse and the longitudinal modes. An explicit choice of the renormalization point will be done later (together 
with the virtual vertex corrections). 

B. General properties of the virtual vertex corrections 

Computing this type of correction is most important for our goals because the vertex strongly depends on the 
renormalization condition and is affected by infrared singularities. It is a very sensitive tool for the physical analysis 
of various renormalization prescriptions and from this calculation we draw our major conclusions. The part of the 
gluon polarization density of the order with the virtual vertex corrections has the following form, 

^2 r j4 7,j4. 



VV 



„/V2 r rl^krl^n 

= 26'^"^ J ^^V"^'^"%q-p,p,-q)VP^fP-{q,k-q,-k)x 

X [-Dret^oi imadvr'^' V^^^^^ {k, -p, p - k)T>*f"'- {k - p)V^-^^' {k - p, q - k, p - q) X 

X [D^r (g)D^r ^ {q - pWiq - k) + D'^^ (9)D?f {q - p)D?f (9 - k)] . (6.19) 
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By virtue of the relations, Dqo = Dg + Di/2 and Du = —Dg + -Di/2, the string, Doo{q)DQo{q — p)Doo{q — k) + 
Dii[q)Dii[q — p)Dii{q— fc), which consists of propagators that form the loop of the virtual vertex, can be rearranged 
into 

Di{q)D,{k - p)Ds{q -k) + D,{q)Di[k ~ v)D,{q - k) + D,{q)D,{k - p)Di{q - k). (6.20) 




FIG. 5. Virtual vertex corrections. The original sum Vooo + Vin can be transformed into Vret + Vadv, thus providing the 
causal behavior of the ladder with the radiative corrections in the order 



In Appendix 1, we show that the vertex function inside ( 3.19| ), rewritten in the form (3.20) is, in fact, the real part 



of the retarded vertex, viz. a triangle graph in which the space-time point x (corresponding to the final momentum 
p) is the latest point on a real time scale. Therefore, one can safely replace the original Vqoo by Vret and Vm (of 
the conjugated graph) by Vadv, since the non-causal remainders do not contribute to the virtual vertex corrections. 
Thus, we proved that the virtual loop corrections do not break the causal picture of the QCD evolution, which has 
been established above, at the tree level. Two retarded virtual vertex loops (depicted at Figs.||c and ||d lead to the 
following analytic expression 

'iN^ r d^k 



n^f-^(p) 



= 6^%^0p'w°,\p)+pyPw°\p)]^S^'>—^ / 2Re^'W^,ir'{k-p,p,-k) 



ewf^\k) w^°'\k) 



'^P2/32 lAii?./; \^^JA^l \ ^P2'^ 



y/32/3/3i(^^_p^p_^P#,/3i^i(^_p) (g21) 



^Wf(fc)Wi^(fc) '''W^{k)Wi{k)_ 

Since the QCD evolution fully develops before the measurement, the physical system of the emitted fields must be 
coherent. These fields represent, at most, a virtual decomposition of the hadron in terms of a set of modes which 
should emerge as propagating fields only after the interaction happens. This understanding of the QCD evolution 
provides us with the natural condition for the renormalization of the fields in the QCD evolution process. Namely, 
radiative corrections to the propagation of the fields should not cause any phase shifts along the line of propagation of 
the hadron. If the hadron is assumed to move along the light cone, the real parts of the vertices and the self-energies 
of the retarded propagators should vanish at p^ = 0. Different subtractions are also explored below. 

The retarded vertex function \ret{p,k ~ p,^k) is proved to be analytic in the lower half-plane of the complex 
light-cone energy, p^ . Equivalently, the advanced function "Vadvip, k — p^—k) is analytic in the upper half-plane of the 
complex light-cone energy, p^ . According to the Schwartz symmetry principle, they form a single analytic function in 
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the entire complex plane with the cuts along the segments of the real axis where the imaginary part of this function 
does not vanish. These segments are pf/p'^ < p^ < oo and —oo < p^ < —{pt — ktY /{k^ — p'^). 

Applying the Cauchy integral formula to the contour depicted at Fig. |^ one can prove the dispersion relation (with 
the subtraction at the point p^ — —51), 

Re^^WiZ\p,-k,k-p) = - du;-{ W3)v(S)(c.,~fc,A:-^) , (6.22) 

where a four-vector w'^ = (jj'^ , lu^ , pt) was introduced for brevity. The subtraction in the dispersion relation (6.22) 
performs an ultraviolet renormalization of the vertex function. The subtracted term is quasi-local and does not affect 
the causal properties of the vertex. The subtraction can only be real and can be performed only at a point where the 
imaginary part is zero. This limits the position of the subtraction point to the segment 

-^P^'-'^'^\p-<4 (6.23) 



fc+ — p+ p~^ 

of the real axis, where Im^'^-'V^et = Im'^^^'Yadv ~ 0. Position = —17 of the subtraction point can be translated into 
the value of p^, p^ = — /i^ = —p^Q — p^. (Since we discuss the vertex function in its natural environment (which 
includes the mass-sell delta-function of the emitted gluon, 6[{k — p)^]) p'^ should be replaced by —k~^{pt — kt)"^ /{k~^ — 
p+).) Under an assumption of strong ordering of the transverse momenta, pt >> kt, the allowed values of //^ are 

0<m'<P?. (6.24) 



© 




FIG. 6. Contour in the plane of complex p used in the proof of the dispersion relations for the vertex function. It is closed 
by two arcs at the infinity. 

This inequality may look abnormal since it does not allow one to choose the renormalization point for the vertex 
function (corresponding to the currently resolved momentum, p), e.g., at p^ — (the much higher momentum which 
is measured at the end point of the QCD evolution) . Such a choice would be almost obvious and even preferable in 
the S-matrix approach to the QCD-evolution which does not require any space-time ordering in the evolution process. 



Nevertheless, in our approach, the upper boundary in the inequality (3.24) is a strict consequence of causality. Indeed, 
the momentum transfer Q will be detected later on, and it cannot affect the earlier loops which correspond to the 
currently resolved momentum, p. Causality forbids one to renormalize what is not yet resolved. However, we are free 
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to choose for the renormahzation any of the aheady resolved lower scales. The limits prescribed by the inequalities 
(I^HI) and are obtained from the dispersive approach which expresses the real parts of the loops (regulating the 

phase shifts) via the imaginary parts (which are related to the real emission processes and are causally ordered in space 
and time). The renormalized real parts of the propagators and vertices will eventually be combined into the effective 
coupling. Thus, the running coupling itself becomes a dynamically-defined quantity that experiences evolution not 
only in terms of the momenta, but also in terms of the space-time variables. This is manifestly in contrast with the 
S-matrix theory which achieves a running regime of the coupling constant by means of the renormahzation group 
method and assigns the same value of coupling to the entire scattering process. 



VII. RETARDED VERTEX IN THE TT-MODE 



The dispersive part of the V^et yields the following term in the evolution, 
"dG(p+,p2)]^^'(^) _ g4^2 r dk+dkt dG{k+,kj) 



TT 



dkf 



p+'S>i{p,k)+p+^2ip,k)] , 



(7.1) 



4(27r)6 J 2fc^ 

where the kernels $i(p, k) and ^2{p, k) come from the dispersion integrals along the right and left cuts, respectively. 

_6+[{k~py] f {p-+n)du;- 
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x8{-2Cip2(^- - p-) - 2C2{k - qf{uj- - p-) - Ci{uj- ~ p-f + C^[{q - kff + C^ip'f ~ 2Cep''{k - qf} , (7.2) 
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^ 8i74(w" -P~f + 8i/5bT + SHebT} , (7.3) 



where the expressions in curly brackets are the result of computing the trace of the product of many Lorentz matrices. 
The functions Cj — Cj{k^ ,p'^ , q'^) and Hj = Hj{k'^^p'^,q'^) are lengthy rational functions. They were obtained 
and further handled with the aid of Mathematica and FeynCalc. The terms proportional to k^ are omitted because 
of the assumption of strong ordering of the transverse momenta along the ladder, k^ « p^. Here we give only an 



explanation of the main steps. As is clearly seen from (7.1) and (7.2), (7.3), the real part of the retarded vertex is not 
calculated and renormalized separately. It is put (in its tensor form) into the environment of the ladder in the form of 
the dispersion integral with one subtraction, and only after that is the trace computed. We had to take this approach, 
because it turned out to be very difficult to find the tensor decomposition of the third rank vertex tensor (similar to 



Eq. ( 3.39 ) for the self-energies). Technically, the calculations are conducted as follows. We start with integrating out 
the variables g^, tj^, and p^, with the aid of three delta functions. Then, we integrate over the angle in the plane 
of the vector gj. All these integrations are finite. The last step is integration over qt — \qt\. Some of the integrals 
happen to be divergent at lower or upper limits of integration, or even at both of them. However, due to miraculous 
relations between different coefficients Cj and iJj, the sums of the integrals are finite. They are as follows: 
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(7.4) 
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(7.5) 

-51. The remaining 



k+{p+-q+)\ k+pl J \ k+{p+-q+)'- k+pi ' 

Recall that the parameter /i^ = is an equivalent of the original subtraction point at p" 

integration over the light-cone component q^ (of the momentum which flows around the loop) is singular at the end 
points. We have carefully checked that no spurious poles connected with the incomplete fixing of the gauge for the 
longitudinal part of the gluon field is involved in the formation of these physical singularities. They all are collinear 
(mass) singularities associated with the end points of the allowed phase space. Moreover, these singularities already 
appear in the calculation of the imaginary part of the vertex function thus being connected to real processes in the 
complementary domain of the external momenta. Hence, we are allowed to introduce a cut-off parameter e (to shield 
the singular behavior), which defines an actual resolution of the process and serves as a definition of a part of the 
newly created gluon field, that should not be treated as real emission. 

Before the integration in Eqs. ( |7.4D and ( |7.5| ), it is expedient to select a large parameter that defines major contri- 
bution to the final answer. The logarithms ln(pf //i^) are good candidates since they can be large if the ratio pi/ is 
either very large or very small. Traditionally, one would wish to associate fi^ with the highest momentum transfer 



in the process. Unfortunately, inequality ( |6.23D deprives us this kind of luxury. Moreover, equations ( |7.4D and ( |7.5| ) 
contain fi^ in the first power, which also makes /i^ = very unlikely. Thus, if we wish to hunt for large logarithms 
of formally, we are forced to choose the renormalization point somewhere near the left tip of the right cut in 

the plane of complex p~ . This corresponds to the on-mass-shell (with respect to the momentum p^ ) renormalization 
of the vertex function, that is, in the domain which is most unphysical in the context of the QCD evolution. We 
even start the perturbative expansion with coupling typical for the non-perturbative domain. Regardless to these 
reservations, we formally obtain pj//^^ 3> 1 as a large parameter and can continue expecting to arrive at the leading 
logarithmic approximation (LLA) with the largest term ~ ln(pf //i'^). The explicit calculation yields 
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(7.6) 



where z = p^ /k'^. The answer occurs to be doubly logarithmic which, may invalidate the perturbative expansion, 
because no terms ^ In^ pi appear in other diagrams which could compensate it in the vertex. 

There is no other way to proceed except for changing the renormalization point. Let us make the new choice 
relying on the causal picture of the process of the QCD evolution, which has been discussed previously, and make a 
subtraction at the point p^ = — 17 = which corresponds to /i^ ~ p1- Then the parameter ^i'^ /pi becomes unity, and 



equations (7.4) and ( |7.5[ ) acqiurc the following form: 
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The remaining integration is straightforward and results in the expression 



(7.8) 
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(7.9) 



It is interesting to trace the origin of the log'^(j>^/e) in this formula. Originally, the calculation leads to three 
doubly-logarithmic terms, 

1 2 ^ 2 k^ , 2 

z m z m h in — , 

e e e 

of which the first term comes from the vertex cut through the momentum p, while the second and the third terms 
come from the cut through the momentum k — p. Even though there is partial cancelation between the first and 
the second terms (which removes two of the double logs), this is not a systematic effect. Indeed, since the kinematic 
regions of two cuts do not intersect, this cancelation cannot be an effect of destructive interference . 

Finally, we have a correction due to the retarded self-energies D^ret'j(p). They are depicted in Fig. ^. The result 
with the arbitrary subtraction point is as follows. 
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(7.10) 



Since the renormalization of the self-energy of the propagator is subject to the same condition as the vertex function, 
we have to take \j? = p\ and thus gain no large Inpj coming from the ultraviolate renormalization. 






FIG. 7. Virtual corrections due to the retarded self-energies of order a^. 
Several remarks are in order: 

1. The dispersive part of the retarded vertex function, renormalized in accordance with the requirement of causality, 
does not contain any -dependent large numbers coming from the UV cut-off. The non-dispersive part of the retarded 
vertex Re '^'■^N^ff^ may undermine the causal picture of the evolution, but it naturally combines with the residue 
ARe ^^^V^^j'', so that in the TT-mode of the evolution equations they together produce the -independent quasi- 
local function, which identically vanishes after renormalization. This is demonstrated by an explicit calculation in 
Appendix 2. The same is true for the terms originating from the four-gluon interaction. They prove to be quasi- local 
as well, and vanish after subtraction at any value of p^ . The proof is also given in Appendix 2. 
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2. The retarded self-energy with momentum p^^ also must be renormalized by the subtraction at the point p~ = 0, 
and thus carries no logarithmic dependence on p^. Therefore, we come to the conclusion, that in the causal picture of 
QCD evolution, the rcnormalization of virtual loops does not produce large logarithms oipf. Hence, the only source 
of these logarithms may be the real processes. 



3. Numerous collinear cut-offs in Eqs. ( |7.9| ) and (7.10) appear already at the level of the calculation of the imaginary 
parts of the retarded vertex and retarded self-energy (unitary cuts) and thus are exclusively due to real processes. 
These singularities are intimately connected with the definition of emission as a physical process and require the 
parameter of resolution. 

4- The expected large logarithms of p^ which are known to drive the QCD evolution in the DGLAP equations 
cannot be an artifact of the rcnormalization prescription. Hence, they have to be outside of the jurisdiction of 
the renormalization-dependent virtual loops and, therefore, must originate from the real-emission processes. The 
rcnormalization prescription itself has to be motivated physically. 



VIII. REAL EMISSION IN THE TT-MODE 



There are four diagrams which contribute to real emission. Three of them originate from the vertex correction (3.7) 
with R ^ S (Fig.|^,b,c). They differ by the place where the initial state density D^^ is built into the diagram. If 
Dqi is put on one of the lines which enters the external vertex of '•^^Hqi we obtain the diagrams which include the 
cut vertex (Fig.||a,b). If D^^ replaces the internal line Dj^s, then we obtain a fragment of the ladder with the crossed 
rungs (Fig. ^jc, the interference term for the emission of two off-spring gluons). The fourth diagram (Fig. ^) is a 
direct descendent of the simple ladder diagram and its order is enhanced due to the further splitting of the decoupled 
gluon. 





C D 



FIG. 8. Real processes in the order a^. 

The two cut vertex diagrams, added together, have a fragment 

Doo{q)Dn{k - p) + Dn{q)Doo{k - p) - ^Di{q)Di{k ^ p) - 2Ds{q)Ds{k - p). (8.1) 

The first singular term on the right side contributes nothing, while both propagators of the second term acquire a 
principal value prescription. The analytic expression for the two cut vertex diagrams is as follows. 
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where Rj — Rj{k^ ,p^ , (7+) are lengthy rational functions. The general line of calculations is approximately the same 
as in the case of the virtual vertex. First, we integrate out the variables q~ and p" using the two delta functions. 
Then, we integrate over the angle in the plane of the vector eft. As previously, these integrations are finite. The 
last step is the integration over qt = \qt\- Some of the integrals are divergent but, due to relations between different 
coefficients Rj the sum of the integrals is finite. It already contains ln(fc+ — (7+); therefore, the final integration dq~^ 
between p+ and fc+ is double-logarithmic. The full answer is, 
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(8.3) 



Here, the term In (p+/e) is a partner of the similar term in the virtual vertex (7.9) which originates from the dispersion 
cut due to the same real process. If a cancelation between these terms had taken place, one could have considered it 
as a result of the interference (which is possible at least owing to the same geometry of the emission.) However, these 
terms not only differ by a factor of two, but also have the same sign. 

The interference diagram (Fig. ^) is calculated in a similar way to the cut-vertex diagram. The initial analytic 
expression is 
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where k = p -\- k — q. The final answer can be cast in a form 
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36z 6z 
Neither the "cut vertex", nor the "interference term" has logarithms Inp^^ but both include logarithms ln(f/.T), which 
are not small at low x. 

The last of the three diagrams related to real processes is due to the decoupled heavy cluster as an emission field in 
the final state. To the lowest possible order, one can mimic it by the diagram in Fig. where the "heavy off-spring 
gluon" decays into two massless gluons. In the lowest order, we have to keep the propagator of the heavy gluon bare, 
which results in the following expression. 
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(fc — pY 

where the cut-off m? , is defined by the properties of the real emission processes, of which the most important is the 
existence of the minimal mass m of the emitted jet, which is of the same order of magnitude as the hadronic mass. 
The two integrations dp'' are straightforward: 
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Thus, we have obtained the first and the only logarithm of p^ which is scaled by a hadronic mass. This scale does not 
coincide with an artificial renormalization scale yu^. This is a real emission process, which produces large logarithms, 
not virtual loops! The final answer becomes 
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Since wc have isolated the mechanism responsible for the large /og's, it is natural to try to elevate its status. This 
does not need special effort, since we have already started from the dressed form of the correlator, which corresponds 
to the emission rung, T)f^^{k — p)nf^(fc — p)^tciy{k — p). This string is nothing but the imaginary part of the full 
propagator {k — p) in the kinematic region A;+ — p+ > 0. Indeed, since the polarization loop 11^^ (g) vanishes at 
positive g+, we can write a chain of transformations, 
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Substituting the last form of this string for the final-state density of the radiated glue, we arrive at 
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To compute the integrals dp~ in this equation we may use the Schwartz symmetry principle. Indeed, basic analytic 
properties of the retarded and advanced gluon self-energies in the plane of complex p~ stimulate similar analytic 
properties of the gluon propagator. Therefore, the real integration can be replaced by the complex integration along 
the path C, which envelopes the cut in the complex plane of the function 



Wj{k-p) = <{ '^iJj' 
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and is closed by a big circle at the infinity. Then the integral is calculated as the residue in the pole of the integrand. 
For example. 
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1 — Rew[^*(fc — p; p~ = Pf/p+)) 

Two properties of the integrand should be noted. First, the pole at the point {k — p)^ = (corresponding to the free 
propagation) is obviously located at the second sheet and does not contribute to the integral (this pole has even to 
be read with the principal value assignment). Second, at the pole p~ = p^/p'^ , the function w{k — p) is real. The 
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other two integrals are computed in the same way, the only difference being, the pole at = is of the second order. 
In the simplest case, when the final-state jet is simulated via a two-gluon system, the functions RewJ'^*(fc — p), which 
enter the final answer are, 
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We have to substitute the value of the argument in the pole, 
ordering, pt^ kt) into Eq, 
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13). Eventually, we express the result in terms of two effective couplings: the transverse 



coupling aT{p1/ z) due to the transverse mode of the emitted heavy gluon, and the longitudinal coupling aL{p1/ z), 
due to the longitudinal mode of the emitted gluon, i. e. 
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The running coupling with the same argument, a{p1/z), was obtained by Dokshitzer et al. by means of the 
"dispersive method" based on an ad hoc dispersion formula for the effective running coupling. Our construction 
(Eqs. (^.9[) ~ (8.15) ) provides a formal proof to this approach. In general, we should not rely on a simple two-body 



model for the density of final states used in these equations. Instead, as is conceived in Ref. [Q, we can consider 
real density of the final state jets which are not yet calculable from a theory, but can be extracted from the data. In 
our picture, with the proton (or nuclei) kept intact before the collision we have to renormalizc the gluon self-energy 
n^(p) at the point p^ = 0. Thus, we have to set p^ = pj in Eqs. (B.14) and (8.15). Then the pt-dependence of the 



effective coupling may only come from the the coUinear cut-off, e, which accounts for the properties of the states of 
emission. 

IX. RADIATIVE CORRECTIONS IN THE LT-MODE 



There are no major differences in the calculations for this mode compared to the TT-mode, since we only change a 
projector. Therefore, we shall omit details in the analysis below. The dispersive part of the virtual vertex before the 
final g^-integration is 
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(9.1) 



The UV-renormalization is performed by means of one subtraction at the point p = — fi. The first and the second 
integrals are obtained via the dispersion relation from the imaginary parts corresponding to the unitary cuts of the 



36 



vertex diagram near the external lines with momenta p and (k — p), respectively. Only the first term of the integrand 
is significant for the final answer. The reason is that the contribution from the cut vertex diagram corresponding to 
the real two-giuon emission, although divergent, exactly coincides, but with an opposite sign, with the second term 
in Eq. ( |9.lD before the renormalization of the latter by means of subtraction in the dispersion relation. Since the 
subtracted counter-term by its definition is quasi-local, the net yield of these two terms will be a quasi-local counter- 
term only. This term is independent of p^ in momentum space and proportional to S{x^) in coordinate space. Hence, 
it must be discarded. With our previous renormalization on the light cone, /i^ = Pt , the answer reads as 
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Unlike the TT-modc, the non-dispersive part of the vertex correction is not quasi-local. After renormalization, it 
looks like 
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With the same choice of the renormalization point at /i^ = p^ we obtain 

Finally, we have a correction due to the retarded self-energies D^rct-j(p). Before specifying the subtraction point, it 
is given by, 
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Since the renormalization of the self-energy of the propagator is subjected to the same condition as the vertex function, 
we have to take /i^ = pi and thus obtain no large Inp^. Thus, as in the case of the TT-mode we can perform the 
UV-renormalization in such a way that no large Inp^ arises. 

Real processes have log^s of pj. One of them, the cut vertex, is gone. The two remaining diagrams correspond to 
the emission of a heavy gluon and the interference between two sequential gluon emissions. 

Let us begin with the heavy gluon. As previously, we start with the lowest order contribution (with the bare 
propagator in the intermediate state). 



. + 2,,.HG ig^N, f dk+dkt dG{k+,kf) 2{2 ~zf _ wf'jk-p) 



Since the density of states 'wf^^{k — p) only establishes the limit of the integration dp^ (being constant above the 
threshold), the integral dp~ requires two cut-offs, the minimal mass to of a gluon jet, and the maximal mass, M, 
resulting in 



r^, + 2^MffG ff^A^e' f dk + dkt dG {k+ , k^) {2 ~ z)\^ . , ^ + ,^ 

[gip+,pm^^=-^ J dfc^ ^(T3^^^(^"-^"'^)l- 



M2 



■7) 



No dependence on p^ appears here, just a big number connected with the end-points of the mass spectrum of the 
final-state jets. This has to be compared with the case of the TT-mode (Eqs. (|8.6D and (|8.7|) ). The difference is easily 
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understood, since the propagator of the longitudinal field has no causal pole. The next logical step is to consider 
the same process with the dressed propagator of the final-state gluon field, as has been done for the case of the TT- 
transition mode. This makes the integral converge; the cut-off M is unnecessary, and we find, 

+ 2^^l^^G ig'^N, f dk+dktdG(k+,k^) 2{2-z)^ 



{k - p)2(l - wl^^ik - p)){l - - p)) 



In our model, the propagator of the final-state time-like gluon has no resonance poles and the integrand of the equation 
(9.8) lacks the propagating pole at p^ = 0. Thus we cannot continue with the contour integration and have to proceed 



in a straightforward way with the simple model of the two-body final state. Then, the integral in the previous equation 
is calculated as follows, 

r^, + 2v,-\HG g^Nc [ dk+dktdG{k+ ,k1) 2{2- zf aictd,rml3T n^ 
[g{p = j ______ . (9.9) 

The pole at z = 1, which was a subject of concern in the order g^, is persistently reproduced in all radiative 
corrections in the order g*, and it requires the physical screening. 



X. SCREENING EFFECTS IN THE QCD EVOLUTION 



We have shown above that the UV-renormalization of the evolution equations can be performed in such a way 
that it does not bring any scale into the problem. Therefore, the only scale which is unambiguously present in the 
equations of the QCD evolution is the one associated with the properties of the final states that may regulate mass 
divergences. In the evolution picture based on null-plane dynamics, this scale systematically enters via the collinear 
cut-off e . Now, we want to find an explicit expression for this cut-off in terms of the parameters of the collective 
modes of the final state. In the simplest case this will be the mass of the transverse plasmon. 

It turns out that the null-plane dynamics, which employs the light-like direction as the Hamiltonian time, 
is ill-suited for computing the screening effects. In null-plane dynamics, the screening is just absent because the 
geometry of fields of these dynamics is singular and cannot provide a mass to a plasmon. In what follows, we prove 
this statement and suggest an alternative type of dynamics, which allows one to derive the evolution equations with 
screening. 

To compute the screening effects, one has to rely on some distribution of the excited modes. One can consistently 
introduce the distribution only in connection with the procedure of its measurement. Throughout this paper, we 
concentrat on the inclusively measured one-particle distribution and argue that the full final state of the collision 
is always prepared as a fluctuation before the moment of the inclusive measurement (and that the measurement 
only excites the system of the allowed final states.) Below, supporting this idea, we construct the basis of states 
that naturally complies with the Lorentz contraction of the incoming objects and includes both strongly and weakly 
localized field modes. In terms of this basis, a snapshot of a virtual fluctuation before a collision (with as of, yet, 
coherently balanced phases) already looks like a distribution of hard particles (with high pt) moving through the 
slowly-varying fields (with low pt). The interaction just freezes this fluctuation with the wrong phases and converts it 
into the true distribution, in which the slow fields are immediately screened by the interaction with the hard particles. 
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Thus, even at a very early time after a collision of two nuclei, the emerging final state is rather a system of plasmons, 
than an ensemble of free partons. 



A. The absence of screening in null-plane dynamics 

The QCD evolution in high-energy processes is always described in the null-plane dynamics, which is the only 
one where the standard variable of evolution equations, the Feynman x, is unambiguously defined. It is tempting 
to try to compute the screening effects staying within framework of these dynamics. The procedure seems to be 



straightforward. One has to compute the retarded propagator (3.41) for the field of emission. 



accounting for the effects of the final-state interaction (collective nature of the final states) in the invariants p^wf (p) 
and w^{p) of the gluon self-energy. In this formula, the first and the second terms are the full propagators of the 
transverse and longitudinal fields of the null-plane dynamics, respectively. The self-energy p'^wf {p) of the transverse 
field has a dispersive part and non-dispersive part. The former is given by the Eq. ( |6.12 ), and it always is zero at 



0. Mathematically, this zero appears in the imaginary part of (p), and survives through the calculation of 



the dispersion integral (6.15). Physically, this happens because the propagation of the field in the dispersive term of 



(T) ( D) 

the self-energy is mediated by the transverse propagator D^J , which maintains the causality of the 11^^/ . This is true 
for all physical gauges, which explicitly separate the transverse and longitudinal fields. 

In the non-dispersive part ( |6.1S ) of Hret (p) , the field propagation is mediated by the longitudinal field, which results 



in Eq. ( |6.1^ ) for p'^w[^^\p). In the case, when the states are populated with density n{k) ~ n(k^, kt) this equation 
is modified to, 

,{ND), _ ~g^N, f ^ 2^_/l + 2n((7) fq++p+Y , l + 2niq-p) f q+ - 2p+ 



This expression has no zero at = and seems to be a good candidate to form the mass term in the dispersion 
equation for the transverse field. However, as before, this term is quasi-local in (it is independent of the light-cone 
energy p^ ) and vanishes after one subtraction at any value of . 

The contact interaction in the tadpole part of the gluon self-energy behaves in the same way, and results in the 

-independent quasi-local term. 

Thus, we can draw two conclusions, (i) The mass of a plasmon with momentum p is generated by the amplitude 
of the forward scattering of a gluon on some distribution of states only, if the gluon field between the two successive 
interactions is longitudinal, or if the interaction is contact, (ii) In the null-plane dynamics, these interactions are local 
in a;+. Therefore, they are incapable of producing any screening effects, which are genuinely non-local. In order to 
incorporate screening effects into the evolution equations, we have to exchange the null-plane Hamiltonian dynamics 
for another dynamics in which the fields would behave in a less singular way. 

B. Field states in the proper-time dynamics 

It was explicitly demonstrated that it is impossible to adequately describe the basic process of forward scattering, 
which is responsible for the generation of the plasmon mass, in the null- plane dynamics. The problem arises due 
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to the singular behavior of the field pattern which is defined as the static field with respect to the time x~^. This 
singular behavior shows that the choice of the Hamiltonian dynamics and proper definition of the field states is highly 
nontrivial and important. The problem we encounter in connection with screening is even more general. Indeed, a 
search for the QGP in heavy-ion collisions is, in the first place, a search for evidence of entropy production. Before 
the collision, the quark and gluon fields are assembled into two coherent wave packets (the nuclei), and therefore, the 
initial entropy equals zero. The coherence is lost, and entropy is created due to the interaction. Though one may 
wish to rely on the invariant formula S = Trplnp, which expresses the entropy S via the density matrix p, at least 
one basis of states should be found explicitly. Thus, it is imperative to find a way to describe quarks and gluons of 
both nuclei as well as the products of their interaction using the same Hamiltonian dynamics. An appropriate choice 
for the gluons is always difficult because the gauge is a global object (as are the Hamiltonian dynamics) and both 
nuclei should be described using the same gauge condition. 

Quantum field theory has a strict definition of dynamics. This notion was introduced by Dirac [p8| at the end of 
the 1940's in connection with his attempt to build a quantum theory of the gravitational field. Every (Hamiltonian) 
dynamics includes its specific definition of the quantum mechanical observables on the (arbitrary) space-like surfaces, 
as well as the means to describe the evolution of the observables from the "earlier" space-like surface to the "later" 
one. 

The primary choice of the degrees of freedom is effective if, even without any interaction, the dynamics of the 
normal modes adequately reflects the main physical features of the phenomenon. The intuitive physical picture clearly 
indicates that the normal modes of the fields participating in the collision of the two nuclei should be compatible 
with their Lorentz contraction. Unlike the incoming plane waves of the standard scattering theory, the nuclei have 
a well-defined shape and the space-time domain of their intersection is also well-defined. Of the ten symmetries of 
the Poincare group, only rotation around the collision z-axis, boost along it, and the translations and boosts in the 
transverse x and y-directions survive. The idea of the collision of two plane sheets immediately leads us to the wedge 
form; the states of quark and gluon fields before and after the collision must be confined within the past and future 
light cones (wedges) with the xy-collision plane as the edge. Therefore, it is profitable to choose, in advance, the set 
of normal modes which have the symmetry of the localized interaction and carry quantum numbers adequate to this 
symmetry. These quantum numbers are the transverse components of momentum and the rapidity of the particle 
(which replaces the component of its momentum). In this ad hoc approach, all the spectral components of the 
nuclear wave functions ought to collapse in the two-dimensional plane of the interaction, even if all the confining 
interactions of the quarks and gluons in the hadrons and the coherence of the hadronic wave functions are neglected. 

In the wedge form of dynamics, the states of free quark and gluon fields are defined (normalized) on the space-like 
hyper-surfaces of the constant proper time r, = <^ — z^. The main idea of this approach is to study the dynamical 
evolution of the interacting fields along the Hamiltonian time r. The gauge of the gluon field is fixed by the condition 
= 0. This simple idea solves several problems. On the one hand, it becomes possible to treat the two different 
light-front dynamics which describe each nucleus of the initial state separately, as two limits of this single dynamics. 
On the other hand, after the collision, this gauge simulates a local (in rapidity) temporal-axial gauge. This feature 
provides a smooth transition to the boost-invariant regime of the created matter expansion (as a first approximation). 
Particularly, addressing the problem of screening, we will be able to compute the plasmon mass in a uniform fashion, 
considering each rapidity interval separately and using the (local) temporal axial gauge. 
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The feature of the states to collapse at the interaction vertex is crucial for understanding the dynamics of the 
collision. A simple optical prototype of the wedge dynamics is the camera obscura (a dark chamber with the pin-hole 
in the wall). Amongst the many possible a priori ways to decompose the incoming light, the camera selects only one. 
Only the spherical harmonics centered at the pin-hole can penetrate inside the camera. The spherical waves reveal 
their angular dependence at some distance from the center and build up the image on the opposite wall. Here, we 
suggest to view the collision of two nuclei as a kind of diffraction of the initial wave functions through the "pin-hole" 
t = 0, z = in <z-plane. 

The states of the wedge dynamics appear to be almost ideally suited for the analysis of the processes that are 
localized at different times t and intervals of rapidity rj and are characterized by a different transverse momentum 
transfer. With respect to any particular process, these states are easily divided into slowly varying fields and localized 
particles. In this way, one may introduce the distribution of particles and study their effect on the dynamics of 
the fields. Thus, we now can calculate the plasmon mass as a local (at some scale) effect which agrees with our 
understanding of its physical origin. 

For the simplest qualitative estimates, it is enough to consider the one-particle wave functions of the scalar field. 
Let us rewrite the wave function ipe.p^ (x) in the following form, 

1 . 0, , ■ ^ , - - f 4-17^-3/2 -imj_Tcosh(j)-e) ipirj_ > f) 

ipe,p^(^) - - 1 4-l7r-3/2g-im_Lrsinh(r,-e)gipj_ri^ < . ^ 

where = m± coshfl, = m± sinh^ (0 being the rapidity of the particle), and, as usual, = p\ . The above 
form implies that r is positive in the future of the wedge vertex and negative in its past. Even though this wave 
function is obviously a plane wave which occupies the whole space, it carries the quantum number 6 (rapidity of the 
particle) instead of the momentum pz ■ A peculiar property of this wave function is that it may be normalized in two 
different ways, either on the hypersurface where t = const, 

[ re'yJx)z^ i,e,pJx)^S{e-9')6ip^-j/^) , (10.3) 

J t— const 

or, equivalently, on the hypersurfaces r = const in the future- and the past-hght wedges of the cohision plane, where 

t2 >0, 

^*e',p'J^) i^9,pA^) - s{e - 0')S{p^ -jf^ . (10.4) 

T— const ^ 

The norm of a particle's wave function always corresponds to the conservation of its charge or probability to find 



it. Since the norm given by Eq. ( 10.4 ) does not depend on r, the particle with a given rapidity 9 (or velocity 
V = tanhS — p^ Ip^\ which is "prepared" on the surface r — const in the past light wedge, cannot flow through the 
light-like wedge boundaries; the particle is predetermined to penetrate in the future light wedge through its vertex. 
The dynamics of the penetration process can be understood in the following way. 

At large m^lrl, the phase of the wave function ipe.pj^ is stationary in a very narrow interval around = 6 (outside 
this interval, the function oscillates with exponentially increasing frequency); the wave function describes a particle 
with rapidity 9 moving along the classical trajectory. However, for to^|t| <C 1, the phase is almost constant along 
the surface t = const. The smaller r is, the more uniformly the domain of stationary phase is stretched out along the 
light cone. A single particle with the wave function ij^e^p^ begins its life as the wave with the given rapidity 6 at large 
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negative r. Later, it becomes spread out over the boundary of the past Hght wedge as r ^ —0. Still being spread, 
it appears on the boundary of the future light wedge. Eventually, it again becomes a wave with rapidity 9 at large 
positive r. The size and location of the interval where the phase of the wave function is stationary plays a central 
role in all subsequent discussions, since it is equivalent to the localization of a particle. Indeed, the overlapping of the 
domains of stationary phases in space and time provides the most effective interaction of the fields. 

The size A77 of the ?7-interval around the particle rapidity 6, where the wave function is stationary, is easily evaluated. 
Extracting the trivial factor e~*™^'^ which defines the evolution of the wave function in the r-direction, we obtain an 



estimate from the exponential of Eq. (10.2), 

2 TOj^rsinh^(Ar//2) - 1 . (10.5) 

The two limiting cases are as follows. 



2 2 

5r] ^ \j , when mj^r 3> 1 , and A77 ^ In , when mj^r <C 1 . (10.6) 

V TOxT mj_T 

In the first case, one may boost this interval into the laboratory reference frame and see that the interval of stationary 
phase is Lorentz contracted (according to the rapidity 9) in z-direction. 



The modes normalized according to (10.4) play the same role as the wave functions of bound states in the examples 
described at the end of Sec. |^ Indeed, the packets of the waves with the same rapidity 9 do not experience dispersion 
in the longitudinal direction. Therefore, an expansion in terms of such waves can be employed to form moving objects. 
The amplitudes and phases of the coefficients in this expansion are balanced in such a way that, before the collision, 
the wave packets represent the finite-sized nuclei with given rapidities. These wave packets must at least partially 
diagonalize the Hamiltonian that includes interactions which maintain the shapes of the nuclei. 

If two localized objects simultaneously pass through the vertex, then the partial waves that form their wave functions 
effectively overlap in the vicinity of the light wedge. The interval of time when the partial wave with transverse 
momentum p±^ is spread out, is of the order r ~ 1/rnx. The high-p^ components of the wave function localize around 
the world line of the initial rapidity t] — 9 earlier, than the \ow-p± components and have less time to interact. If no 
interaction occurs at sufficiently small r, then the amplitudes and the phases of these waves remain unchanged and 
the packet assembles into the initial finite-sized object. 

If the interaction takes place and the balance of phases becomes broken, then the estimates of Eqs. (10. f; ) and ( 10. 6| ) 



show how to construct the picture of screening at sufficiently small r, even before the secondary collisions come into 
the game. The states with pt ^ vary only slowly at the rapidity interval At] ~ ln(2/pfT) 3> 1, and cannot be 
considered to be particles. At the same time, the states with kt ^ r^^ are well localized in the rapidity direction. 



For these states, Srj ^ ^J^jktT <C 1 and we may safely view them as the particles with the distribution which can be 
measured inclusively at the time r. The forward scattering of the particles with transverse momentum kt on the field 
with the momentum pt <C fcf results in an adjoint mass of the soft field which screens the possible collinear singularity 
in the evolution equations. The process of evolution naturally saturates at the proper time r when the adjoint mass 
becomes comparable to the transverse momentum pt ~ 1/r. This is the end of the "earliest stage". At later times, 
collisions between the partons-plasmons take over. 

Another important question is how early the "earliest stage" begins. It is easy to understand that the shortest time 
scale that may be seriously considered in the theory is defined by the full energy of collision, Tmin ^ '^1 \fs- (For nuclei 
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collisions, Tmin where s is the invariant mass per participating nucleon.) At this time, the stationary 

phase of a partial wave is stretched over the widest rapidity interval Srj ^ ln(^/s/m^). It is not surprising that this 
estimate coincides with the well known kinematically allowed width 2Y of the rapidity plateau, 2Y sa log{s/m?^^^^). 
We just have two complementary ways to obtain the same quantity. 

The entire design of pQCD is aimed at getting rid of logarithms like log{s/m?)\ they lead to mass singularities. Their 
presence invalidates the power counting in the calculation of the scattering amplitudes and makes the dimensional 
analysis (including the RG methods) impossible. On the contrary, in AA-collisions, this logarithms play a major role 
defining the energy density in the emerging dense system and, consequently, the screening effects. However, even in 
this case, we may find subprocesses (with the largest pt-transfer), which are insensitive to the collective effects, and 
where the two strategies (screening and removing) must be not too far away from each other. 

A full analysis of the QCD evolution and the collective effects in the scope of the wedge dynamics is a subject of 
separate study which is now underway. To treat these problems in a global fashion, one has to derive propagators 
of perturbation theory in curvilinear coordinates, separate longitudinal and transverse gauge fields, and perform 
quantization of all fields Here we shall address only the estimate of two most important effects: 

(i) The plasmon mass, which may be computed with reference to the local properties of the wedge dynamics. 

(ii) Screening in the evolution equations at large pt, where the difference between the wedge and the null-plane 
dynamics is expected to be minimal. 



C. Local screening and the mass of a plasmon 



The local parameterization of the dynamics based on the system of the space-like surfaces r = const is easily 
implemented with the aid of the vector 1/^(7/) = (cosh 77, Ot, sinh?]), which is normal to the hyper-surface and indicates 
the local "time" direction. As long as we are interested in the process, which takes place in a limited slice in rapidity, 
we may consider 77 as a label of this slice and perform all calculations in the Lorentz frame which moves with the 
rapidity -q. In this frame, u'i — (1,0,0,0), and the global gauge condition, u{rj)A = A"^ = , becomes just A^^ — 0. 
A dynamical formation of the plasmon mass is exactly the case when we can proceed in this simplified manner. 

The bare propagator in this (local) temporal-axial gauge is of the following form, 

1 / k'^v'' k'^ii'^ k'^k'^\ d^^\{k) 

D'-'iu, fc) = -i - g^"^ + " +f -iA^]= + , (10.7) 



k^\ ^ (ku) {kuf J fc2 {ku) 

where the first and the second terms in this equation are the propagators of the transverse and longitudinal fields, 
respectively. The polarization matrices are, 

d^,^,{k)^-g^'^ + u^u^^dll^{k), d^,l,{k)=''^^—^^^^^^^P^, (10.8) 

and the gluon self-energy H'"' is of the following form, 

H^ (p) = d![rj.^(j>) wt{p) + p2 [(p^)2 _ p2] wl{p). (10.9) 

Since H^'' always appears in an assembly DUD the terms, like p^u'' -I- tt^'p" or u'^u'^ , which are necessary to provide 
transversality of H'^", cancel out. The invariants wt and wj^ can be found from the two contractions, 

- <?^,(p)n^''(p) = 2 wt{p) + WLip), and w^u,(p)H''''(p) = ^P^^^^^^w^p) ■ (10.10) 
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The solution of the Schwinger-Dyson equation for the retarded propagator can be cast in the form 

which exhibits an exphcit separation of the transverse and longitudinal propagators. 

Our minimal goal is to estimate the mass of the transverse plasmon (the collective mode) which is expected to 
replace the final state of free propagation in the QCD evolution equations. Therefore, we need to compute the 
self-energy Wj^{p) of the soft transverse mode accounting for the effect of the hard part of the (already measured) 
final-state inclusive gluon distribution, dNg/dkfdy. To pick up the dominant effect, we have to compute w^{p) at the 
point p2 ^ |]. Therefore, we can parameterize the momentum p of the soft mode as p'^ = (pt cosh 9, pt,Ptsmh 9). 
The self-energy w^{p) has dispersive and non-dispersive parts, which are defined by equations ( 6.12| ) and ( |6.13 ), 
respectively. Even before the loop integral is computed, it is straightforward to check that the dispersive part is 
proportional to p^ and can be disregarded in our simple estimate. The non-dispersive part of wi^{p) is represented by 
the integral, 

) g^N, pI f d^q (S{q')..„ , 

^ ^ y ^^^^ 

(go +po)^(5g§ - 2goPo +Po) 9o+Po 



(p- q 



2 



4 

(q^p-q)} + 0{p'), (10.12) 



4(p-q)2 2 

where the "normalization volume", T^^^-'(t), for the hard modes with transverse momenta qt S> pt is defined as the 
volume of the domain on the hyper-surface r = const, where the wave function of the soft plasmon mode is stationary, 
and is given by 

F^^) (r) = 7ri?2 T A77 , A7/-ln— , 

PtT 

where R is the radius of the colliding nuclei, -kR^ w ■nr'^A^/'^, and rp is the proton radius. The expression ( 10.12| ) 
includes a vacuum part which is UV-divergent and requires renormalization (which we neglect for now). Integrating 
go out and setting = we arrive at 



(ND), . _ g^Nc f d\ d N, 



27ri?2 T Ar] J |q|3 d^ q 



(p-q)^ ^ (5|q|^ + 2|qp|pp + |pn ' 

i2 q p 



(10.13) 



2 2(p - q)2 

As for the distribution d Ng/d^q — d Ng/qtd^qtdy , we shall assume that it is homogeneous in the rapidity interval 
A77 centered at the rapidity 9 of the soft mode, n = n{qt,9 ~ y) fa n(qt,9). Furthermore, we can take 9^0 and 
consequently, « 0, po ~ Pt- With the same accuracy, we may take q^ — qt sinhy w 0, dq^ = qt coshy dy « qt dy 
and integrate over qt with the condition qt ^ pt- In this approximation, we obtain 

^r'lP;,";/^")^^ r^im^- (10-") 



7ri?2 J ^2 rqt d qf 



where the function dNg/ dq^ is the density of the final-state gluons in the unit interval near rapidity 9. In the same 
way we compute the contribution of the tadpole term, 

.4«(rf..;,= .0,.^ r^(i,M^)^'-S^r^i^ . (10-15, 



ttR^t Arj J \q\ \2 J ttR^ Jp2 rqt d qj 

where the divergent vacuum term vanishes after UV-renormalization 
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Now, the renormalization can be performed in the usual way by subtracting the quasilocal terms after the integration 
with a weight function n{6). Indeed, the point of renormaUzation can be chosen in the domain 9 =poo {p^ — > 
or ) where it does not affect the finite part of w{p)t produced by the physical distribution n{6,pt) which 

vanishes for modes with infinite rapidity. Assembling these two contributions, we obtain, 

m%{plO) = wr\plO;p- - 0)+4'''\plO-y -0)^(3 + 1)'-^ r Ml^m^ . (10.16) 

■""-ft Jp2 TQt a 

This quantity is defined locally in rapidity and it works as a feed-back that limits the possible energy of the emission 
field from below and thus screens the mass singularity in the evolution equations. 

Because the product ptT in Eq.( 10.l"6| ) is close to unity, the proper time r and momentum pt are the complementary 



parameters, this equation can be read in two ways. On the one hand, the smaller pt is, the more particles with higher 
momentum interacting with the soft mode are emitted, and the heavier this soft mode is. On the other hand, we 
know that the mode with momentum pt acquires the status of the emission field only by a time r ~ The later 

this happens, the more that particles which participate in the formation of the plasmon mass are created, and the 
larger the mass of plasmon is. 

XI. EVOLUTION EQUATIONS WITH FINAL STATE SCREENING 
A. Screening in the order a^. 

Our last problem is to show that the replacement of the free radiation field by the collective modes of the true final 
state screens the collinear singularities that we encountered in the null-plane dynamics. Previously, these singularities 
were coming either from the poles of the polarization sum d^", or from the integration dp^ of the (5+[(fc — p)^]- In 
order to smooth out this behavior, we consider the evolution equations in the same dynamics which allows for the 
effect of local mass generation. While the calculation of the plasmon mass could be performed in the co-moving 
reference frame, the evolution must be treated in a more global fashion. In order to find the space-time domain 
where the QCD evolution described by the DGLAP equation develops, let us go back to the general four-vector 
u'^ = (cosh?]. Of, sinhry), that defines the local time direction at the rapidity rj and keep the latter as a parameter. 
The gauge of the gluon field is fixed by the condition uA ~ which may be conveniently rewritten as 

u{7j)A = ^e-''A+ + ^e'^A- = . (11.1) 

From this, we infer that the gauge condition A^ = holds as i] — s- ~oo, that is, in the nearest vicinity of the 
hyper-plane = 0. In the process of e-p DIS, the Lorentz frame moving with this rapidity is the rest frame for 
the electron, and the infinite-momentum frame for the proton. The momentum of the virtual photon of the DIS has 
components g"*" = and q~ = 2i//P+. Thus, the DGLAP equations describe the process which is localized in the 
space-time domain of the electron fragmentation. 

In what follows, we attempt to answer only two questions: 

(i) Does the plasmon mass shield the singularity in the evolution equations? 

(ii) Does the equation with screening allow for a smooth transition to the DGLAP equation when the screening 
mass becomes small? 
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For these limited goals, we can drop all regular terms. To derive the evolution equation, we start with Eq. (|5.3|) 



and extract the scalar equations by means of Eqs. (10. IC). In this way we obtain an analog of the Eqs. ( |5.4| ), ( |5.8| ) 
for the evolution of the transverse fields, 



dG{x,pf, 



-Nra, 



dpi 



dpo 



2uj 



-d{ko ~po ~uj)- 



2\ ^2 r 

P (—\ 

^99 \ 7,n I 



1 



p" 



fcO [p- 



,212 



where only the terms which do not vanish at s oo (rj — oo) are retained . In this equation, 

= {K - p,)2 +m2((fc; -pl)2) = {k'^~plf + ih -p-f+Tn^ih-p-tf) , 
and the temporal and longitudinal components of momenta are those of the infinite momentum frame, 



(11.2) 
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(11.3) 



Both are large and scale with the momentum of the proton. (In Eq. (11.2), and in what follows the label * is dropped.) 
The operational definition of the structure function in the new variables changes slightly with respect to the previous 



definition (5.7) 



dpQ 



Po w'{!^\p) 



(p2 _ w^{p)){p'^ - wq,{p)) 



dG{x,p1) 
dpi 



(11.4) 

The mass 



where by virtue of ( 11.3 ), the domain of the integration dp^ is strongly localized near the point p° p 
m?({kt — pt)^) was computed in Sec. XC and is used as the pole mass of the emitted gluon-plasmon. It has the 



property to increase with decreasing transverse momentum. Integrating out p^ with the aid of the delta-function in 
the Eq. ( |ll.2D , we arrive at 



dG(p\pf, 



dk,d^ktdG{k',kf) {k^-ujf 



2uj 



dkl 



p2[(fc,-C^)2-p2] 



kz 



uj(kz - ^) 



h.2 



(11.5) 



where the expression in the square brackets is just the splitting kernel Pgg and the strong ordering (p^ ^ fc2-j jg 
assumed from now on. Once again, this equation is written in the first approximation with respect to the small 
parameter pl/k^ ~ Pt/s, and the former pole at p° = k'^ is translated into the non-singular factor l/uj. 

By examination, Eq. ( |11.5 ) has no singularity in the remaining integration dkz . It is reliably screened by the finite p^ 
regardless of the presence of the plasmon mass. Such a screening happens because now the proton has finite rapidity] 
and the dynamics are free from the artificial singularity of the null-plane. However, if we consider a non-ordered 
emission, or just take pt to be small, then only the mass, m{{kt — pt)^), would provide the desired protection. Only in 
the formal limit m? / s 0, the singularity at k^ = Pz becomes a reality, and we recover the usual DGLAP equation. 
It is easy to demonstrate, that the same mechanism of shielding of the mass singularity works for Eq. (|]^), which 
describes the transition between the longitudinal and transverse modes of the gluon field. 



With the mass computed according to ( 10.14 ), Eq. ( 11.5 ) can be considered as the simplest prototype of the 
evolution equation with screening for the fluctuations in the AA-collisions. One should keep in mind, that (by virtue 
of (11.3)) the set of limits, kz — Pz 0, fc+ — p+ 0, and t] —oo, which has been employed for its derivation 
is intrinsically ambiguous. In order to obtain the evolution equations that can handle the emission in the central 
rapidity region we cannot approach the limit of rj s- — oo. These equations can be derived only in the framework of 
the "wedge dynamics" [|o|j2|]. 
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B. Suppression of the radiative corrections in the order q^. 



The last thing that we are going to demonstrate with our Umited "local" approximation is that the radiative 
corrections of the order are doubly-suppressed by the plasmon masses. Indeed, a review of our previous calculations 
in this order shows that all of them are connected with processes with at least two real emissions. An example of 
such a process with two consecutive emissions is depicted at Fig. ^. 




FIG. 9. Fragment of the ladder with two rungs of emission. 
In any of these cases, we have to integrate the expressions with two mass-shell delta functions, 

d'^k S+[ik-p)^ -ml] S+[il- k)^ -ml]Pi{p+/k+)P2{k+/l+) , (11.6) 



where Pi^2{z) may have a pole at z = 1, mi = mi{{kt — Pt)^), = 'm2{{lt ~ ^t)^); ^'^'^ the total emitted momentum 
IS I — p — s. For example, let the coUinear singularity emerge when the integration dk'^ reaches one or both ends of 
the interval p^ < k^ < . Introducing the new variables, f = k — p, we can integrate out f" with the aid of one of 
the delta- functions and remove the second one by integrating over the angle between the transverse vectors ft and si- 
The latter integration defines the limits for the magnitude of the transverse component /i, which has to be confined 
between the two positive roots of the equation, 

list + ft? - (.+ - /+)(.- - ^-^j^) - ml] [(.+ - /+)(.- - ^-Ij^) -ml- ist - ft)'] = . (11.7) 

These roots exist (the measure of integration dft does not vanish) only if the variable f~^ is within the following limits. 



m? — mo / m? -I- mS /m? — m,n \ ^ „ f ^ 



m? — mn / m? + m, / m? — m^ , 

This inequality shows that the domain of integration over fc+ gets smaller with the growth of the masses, and it just 
vanishes when, e.g., either mf > or ml > s^, and the emission of two gluons becomes impossible. Hence, the finite 
plasmon masses not only screen the mass singularities, but also strongly reduce the phase-space for the terms of the 
order a?. 



47 



ACKNOWLEDGMENTS 

We are grateful to Yu. Dokshitzer, A. Kovner, L. McLerran, B. MuUer, and H. Weigert for many stimulating 
discussions. We wish to thank Edward Shuryak for helpful conversations at various stages in the development of this 
work. We are grateful to Jianwei Qiu for attracting our attention to his paper We appreciate the help of Scott 
Payson who critically read the manuscript. 

This work was supported by the U.S. Department of Energy under Contract No. DE-FG02-94ER40831. 



XII. APPENDIX 1. CAUSAL PROPERTIES OF THE EVOLUTION LADDER 

In the course of computing the (7^-elementary ladder cell we met various types of radiative corrections with retarded 
order of their temporal arguments. This circumstance allowed us to make a firm conclusion (which was verified at the 
tree-level) about the time ordering in process the QCD evolution process. It is also fortunate in one more respect. The 
retarded (advanced) functions have the unique analytic properties that they are regular in the upper (lower) half-plane 
of the complex energy. Since we use Hamiltonian dynamics with the coordinate chosen as the "time" direction, 
the corresponding conjugate variable is the light-cone energy, , and we can rely on the analytic properties in the 
complex plane of p^. Beyond the tree-level calculations, we include radiative corrections with loops, self-energies and 
vertex functions. We must check whether they maintain the aforementioned causal behavior. Self-energy corrections 
to the propagators are of an explicitly retarded type, and they do not cause any problems. The loops of the vertices 
have yet to be examined. The two vertex functions, Vboo: f^nd Vm , which come from the terms with i? = 5 in 
Eq.(|6.19D, have the same surroundings and appear in the combination Vboo + Vm. This sum forms the following 



string of propagators which constitute the loop, 

Doo{q)Doo{k - p)Doo{q ~k)+ Dn{q)Dn{k - p)Dn{q - k). (Al.l) 
Since Dqq = Dg + Di/2 and Dn = —Dg + Di/2, this string can be identically rewritten as, 

Di{q)Dsik - p)D,{q -k)+ Ds{q)Di{k - p)Ds{q - k) + D ,{q)D s{k ~ p)Di{q - k), (Al.2) 
where we omitted the term DiDiDi/2, which identically vanishes in the surroundings of the evolution ladder. We 



want to show that expression ( Al.2 ) is nothing but the real part of the retarded vertex, which is an analytic function 
in the lower half-plane of energy p^ . To prove this statement, we should derive the formula for the retarded vertex 
to order ff^ . This is more easily done in the coordinate form. 

Let us start with the definition of the retarded vertex in the form of the functional derivative of the retarded 
self-energy with respect to the field A{z), 

^^^UJ^.y.^ = - ''^'Sy:^ ■ (A1.3) 

g5Aij{zB) 

The self-energy liret is built from the propagators Dret^ Dadv, and Di. The simplest way to differentiate them is to 
use the equations ( |3.5[ ) for Dret, Dadv, and Di in the presence of the classical external field. These equations can be 
written as follows: 
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T>iix,y)^Di{x,y)+ J dz[Dret{x, z)V{z)A{z)Di{z,y) + Di{x, z)V{z)Aiz)Dadviz,y) . 
It is straightforward to obtain, 

^^WZ^{x,y,z) = -'-{[ ^,^)^"^""^ Wi^.TM^,2/)^''^''''K2/)^r"^(2/,^) 

+ ^a,aa. (^)^^".P. (^^ 2^)V^P.PPi (y)15Pi-2 (^^ ^^ya.aa, (^)i^-i^i (y, 3;) 



(A1.5) 



(A1.6) 



Here, the groups of the first three, and of the last three terms differ only by the direction in which the arguments are 
going around the loop. Therefore, these two groups are identical. In the momentum representation, we obtain 

(3)v„t(p, -fc, k-p)^-'-J dxdydze'P^^~'^+'''^y-''^ [ V{x) Dret{x - y)V{y)Dret{y - z) V{z)D,iz - x) 
+V{x)Dret{x - y)V{y)Di{y ~ z)V{z)Dadv{z - x) + V{x)Di(x - y)V{y)Dadv{y ~ z)V{z)Dadv{z - x) ] . (Al.7) 




z + 





FIG. 10. Three diagrams contributing to the retarded vertex. Lines with arrows correspond to the propagators Dret and 
point X is the latest. The crossed lines correspond to the correlators Di, the densities of real states. 



The elements which are responsible for the analytic properties of the retarded vertex with respect to momentum p 
are underlined in the formula and marked by the time-directed arrows in Fig. hol In the same way, we may introduce 



an advanced vertex function and put it in a form similar to (Al.7) 



^""^Vadvip, ~k, k-p) = --j dxdydze'P^^'^-^+^'^^y-'-^ [ V{x) Dadv{x - y)V{y)Dadv{y ~ z) V{z)Di{z - x) 
+V{x)Dadv{x - y)V{y)Di[y ~ z)V{z) Dret{z ^ x) + V[x)D^{x - y)V{y)Dret[y - z)V{z) Dret{z ~ x) ] , (A1.8) 

which shows that at the real light-cone energy p~, Vadv is the complex conjugate of Vret- This is exactly the sum 
Vret + Vadv = Vqoo + Vm (twicc the real part of Vret) that defines the vertex loop correction to the elementary cell of 
the ladder. Indeed, since 



Dret = Ds + 2^07 Dadv = Dg — l^^^i 



(A1.9) 



we find that 
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Re(3)v(2^(p,-fc,fc-p) = -^ J ^[vMq)VyDs{q-k)vMq-p) 

+VMq)VyDi{q^k)V,D,{q~p) + V^D,{q)VyD,{q~k)VMq-p)] . (ALIO) 



Thus, it is natural to use the dispersion relation (3.22) to compute the real part of Vret- However, one should 
remember that each of the retarded propagators of the gauge field includes a part that corresponds to longitudinal 
fields, which are not driven by the Hamiltonian equations of motion. Only the transverse fields propagate and obey 
causality. Therefore, those retarded and advanced functions, which maintain the causal properties of the vertex (and 



are underlined in Eq.(A1.7) ) have to be split into longitudinal and transverse parts, 

Dretiq) = D'^Jiq) + D(^Hq), ^ oZliq) + D^^\q) , (Al.ll) 

and only the causal transverse part has the proper analytic behavior, which results in dispersion relations linking the 
real and imaginary parts of the vertex function. Thus we encounter a necessity of separating the dispersive part of 
the vertex, V''^\ from the non-dispersive part, 

For our practical goals, we need only the real part of the dispersive term of the vertex function, which is, 

(3)v(S'(p,-fc,A:-p)^-^ j ^^[V^D^Zli<l)VyDZl{q-k)VMq^p) + 
+V^Dret{q)VyDiiq - k)V,DZl{q - p) + VM^WyDadviq - k)V,DZl{q-p) ]■ (Al.12) 

It can be computed by means of dispersion relation via the imaginary part, which consists of two terms. The main 
term, is symmetric with respect to the singular functions (commutators Dq and densities of states Di), and is given 
by, 

Im(3)vl-)(.,-fc,.-.) = 4/^x 

X [ VMl)VyDf\q - k)VMq -u:)- VMq)VyDf\q - k)VMq - a>) ] + 
+ [ V,DP{q)VyDo{q - k)VMq ^ u) ^ V^D^ {q)VyDi{q - k)VMq - to) ] + 

+ [ VMct)VyDi{q - k)V,Df\q -uj)- VMlWyDoiq - k)V,Df\q - Cc-) ] j . (A1.13) 

Symmctrization has been achieved by an extra splitting of the two retarded propagators which do not depend on the 
momentum p and do not participate in maintaining the analytic properties with respect to the (light-cone) energy 
p^ . Therefore, we have an additional residue of the splitting procedure which contains the longitudinal propagators 
D^^\ i.e., 



Aim (3)v(f,)(a;,-fc,fc-c.)^-l J- 



(27r)4 

[ - V.Di'^HqWyDiiq - k)VMq - ^) - VMdWyDi'^Xq - k)VMQ - ^)] • (A1.14) 



In the last two equations, wc have introduced the four-vector — [p'^ , uj^ , pi) . The first two terms in Eq. ( Al.13 ) 
correspond to the unitary cut of the diagram for Vret near the external line with momentum p. The next two terms 
correspond to the unitary cut near the external line with momentum p — k. At real w^, the imaginary part due 
to these terms differs from zero for lu^ > Pt/p'^ and for uj^ < —{pt ~ kt)'^ /(k^ — p^), respectively. Therefore, the 
subtraction point in the dispersion relations with respect to the light-cone energy p^ is confined to the finite interval 
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of p~ between the tips of the two cuts. The last two terms in Eq.(A_1.13) would have corresponded to the unitary cut 
near the external line with momentum fc, if the momentum k were time-like. As long as dispersion relations do not 
involve the momentum k, these terms vanislie in the ladder kinematics, since we always have k^ < 0. 

The non-dispersive part of the vertex is not causal and can only be calculated in a straightforward manner, 

(^)viT{p,-k,k-p) = -'- I ^^[V^Di^\q)VyD^^Hq~k)VMq-p) + 
+V^D^''\q)VyD^P{q~k)VMq-p) + V,D,{q)VyDi{q~k)V,D'^'^\q-p) + 

+VM<l)yyDs{q - k)V,D<^'^\q ~ p)\ . (Al.15) 

We use the vertex function only when it is embedded in the ladder diagram, and therefore, we do not need to compute 
it separately. 

XIII. APPENDIX 2. THE VANISHING CONTRIBUTIONS 

There are several types of terms that do not contribute to the evolution equations for various reasons. 

A. The terms with the four-gluon vertices 



In the order g"*, we have eleven diagrams with four-gluon vertices. They can appear only in the TT-transition mode. 
Ten of them turn out to be zero. These are diagrams (a-e) in Fig. |ll| and their "mirror reflections". 








FIG. 11. Diagrams with four-gluon vertices. 



Diagrams a-c include the virtual loops that are the real parts of the retarded functions. Hence, they can be 
computed via the dispersion relations. Even though the dispersion integral is found to be UV-finite, it still requires 
renormalization. Therefore, the dispersion integral has to be rewritten with the subtraction at some point p^ — — fi. 
For the diagram (a) we get, 

"dG'(]5+,p?)] _ jg^^c f dk+dkt dG{k+,kl) ^^^^ f ff duj- duj' 



dp, 



(a) 



{2n) 



2k+ 



dkf 



dq+d''qt5+[{uj-qf]5+[q^]({uj--p-) [k - q) 



[p2 



dp 
p^{2k+q 



uj — p 



q+p+) 



p^ 



k+{p+-q+) 



k+p+2(p+ -q+) 



= 0. 



(A2.1) 



This virtual correction identically vanishes after subtraction at any value of the light-cone energy. Therefore, in the 
coordinate space this term is quasi- local ( ~ S{x^) ) and does not contribute to the evolution. 
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Diagram (b) behaves similarly: it is finite and we calculate it using the dispersion relation over p making a 
subtraction at some point = — il to renormalize it, 



dp't 



{b) 



duj 



ig^N, f dk+dkt dG{k+,kj) . 2. _ 

J2^J^ dk^^^^^-P^^'^P 

diu- \ f dq+d^qt5+[{q-Uj)^]S+[{k~q)^] 



LU — p 



{UJ -p ) 



2 



[p2 



V 



0. 



(A2.2) 



\k+ - q^)(k+ - p+) k+{q+ -p+){k+ -q+) {k+ - p+)k+ {q+ - p+) 
Without renormalization, these two terms would produce finite corrections comparable to those coming from the 
three-gluon vertex (they include both the Coulomb logarithm, and the coUinear cut-off). 
Diagram (c) is oi the form, 



dG(p+,p2) 



dp 



ig^Nc f dk+dkt dG(k+,k^) f,^6\{k-p)^] f dq+d'^qt 
' 'dp ' 



(c) 



2k+ 



dkf 



[p2 



6[{k~qr] d[q'] 



(fc - qy 



{k-qf{LO-p-) [{k-qff p^{k-qf{2k+q+-k+p+-q+p+) 



0. (A2.3) 



p+ k+{p+-q+) k+p+'^{p+ - q+) 

It proves to be zero after delta-functions are used to integrate out p" and q^ . (The expression in the curly brackets 
is zero.) 

The diagrams (d) and (e) correspond to real processes. They vanish identically, which is seen only after the 
integrations dp^, dq^ , and dqt are carried out. 

The only term which has a four-gluon vertex that survives, is given by diagram (f). It corresponds to the real 
emission 



dG{p+,p]) 
dpi 



"l2^ 



(/) 



dk+dkt dG(k+,kl)^^_ dq+d^qt S [jq ~ p)^]S [jk - qf] 

—2 



2fc^ 



dkf 



-dp 



g^N^ [dk+dkt dG{k+,k^) 6(fc+-p+) 



{2ny 



2k+ 



dkj 



k+p+p1 



(A2.4) 



B. The non-dispersive and relative terms in three gluon vertex. 



The two terms in the three-gluon vertex in the TT-transition mode, which were not discussed in the body of the 



paper, are given by Eqs. ( Al.14 ) and ( Al.15 ). To compute the contribution of the longitudinal (non-causal) fields in 
the dispersive part we have to substitute ( A1.14| ) into Eq. ( |6.22| ) and then proceed with Eq. ( 6.21 ). These steps result 



dp\ 



4(27r)6 

where and are given by 

duj^ 



2fc^ 



dkf 



ALl 



duj 



dq+dq- + + r 5[{q-Ujf]5[q^ 

sign((7^ - p^) I d qt — 



LU — p U) + J j 2 " • ' ' ' J Jp2 

{uj--p-) {k-qf p^{2k+q+ ~ k+p+ - q+p+) 



p^ 



k+{p+ -q+) 



k+p+^{p+-q+) 



*l(q^ 



(A2.6) 
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and 



L2 



duj 



dq^dq 



sign((7+ -p+) / d qt 



,5[{q-Lof]5[{k-qf 



(A2.7) 



ik+ - q+){k+ - p+) k+{q+ -p+){k+ -q+) {k+ - p+)k+ {q+ - p+) 

where 5*1, 2(9^) are rational function. After the integration over q^ , uj^ , and both <i>ALi and <&al2 vanish for 
an arbitrary subtraction point fJ. Therefore, the longitudinal fields contribute only a quasi-local term to the vertex 



function. The non-dispersive part, given by Eq. (A1.15), being substituted into Eq. (6.21), results in 
dG(p+,p2)]^^'(^^) g^Nl fdk+dkt dGik+,kf) 



dp: 



TT 



4(27r)6 



2fc^ 



dk^ 



with ^NDi and ^nd2 given by 



9ndi — a q - 



[p 



212 



dp-[p+'S^NDiip,k)+p+^ND2{p,k)]6[ik-pf] , (A2., 



(A2.9) 



p\2k+q+ - k+p+ - q+p+) - (fc - qrip+Y) 



and 



ND2 



= A 



5[{q-pr]S[{k~qy 
[P'Y 



p'ik+-q+)~q'ik+-p+)}^2{q^) 



(A2.10) 



As with all other constituents of the vertex function, the non-dispersive part, even being finite, requires renormal- 
ization. We make a subtraction at p^ = — f2. This results in both functions, ^nd2 and ^nd2, separately being 
zero. Thus, all terms associated with the longitudinal fields in the vertex function in the TT-transition mode are 
quasi-local. 
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